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dynamically  relevant  conservation  laws  governing  failure  and  energy  release/absorption  at  the  failure 
boundary  and  in  the  surrounding  medium,  so  that  in  a  fundamental  sense  these  models  are  ad-hoc. 
Further,  examination  of  theoretical  representations  of  the  seismic  radiation  from  earthquakes  that  employ 
the  divergence  of  changes  in  the  stress  within  the  medium  as  a  source  term  in  the  equations  of  motion 
(such  as  the  "stress  glut  model"),  are  shown  to  be  without  a  rational  physical  basis  and  to  have  no  mean¬ 
ingful  relationship  to  the  dynamics  or  kinematics  of  spontaneous  failure. 

In  the  present  study  we" develop  an  approach  that  incorporates  the  (nonlinear)  conservation  relations 
on  the  failure  surface,  as  well  as  those  appropriate  in  the  surrounding  linear  zone,  to  generate  a  Greens 
function  integral  equation  describing  both  the  failure  growth  and  the  (interacting)  seismic  radiation  field. 
The  method  involves  the  explicit  decomposition  of  total  stress-displacement  fields  into  dynamic  and 
equilibrium  parts,  with  the  latter  dependent  on  time  because  of  the  growth  of  a  new  boundary  correspond¬ 
ing  to  the  failure  zone  boundary  within  the  prestressed  medium;  with  this  boundary  growth  necessitating 
a  time  dependent  readjustment  of  the  prestress  state  to  maintain  equilibrium.  It  is  shown  that  the  time 
changes  in  the  equilibrium  fields  outside  the  failure  zone  give  rise  to  an  "equivalent  force  term”  in  the 
equations  of  motion  in  the  linear  zone  corresponding  to  the  inertial  effect  of  time  dependent  changes  in 
the  equilibrium  displacement  field.  Thus  the  "proper"  equivalent  force  term  in  the  equations  of  motion  in 
the  linear  zone  is  this  inertial  term  and  not  the  divergence  of  the  stress  field,  either  outside  or  inside  the 
failure  zone.  The  Green’s  integral  equation,  arising  from  the  equations  of  motion  and  boundary  condi¬ 
tions  (conservation  relations)  on  the  growing  failure  surface,  show  that  there  are  three  source  effects  that 
mitigate  the  seismic  radiation  field,  in  particular  a  primary  term  resulting  from  the  inertial  effects  involv¬ 
ing  the  equilibrium  displacement  in  the  medium  surrounding  the  failure  zone  (a  volume  relaxation  effect) 
plus  secondary  scattering  effects  from  the  failure  surface  and  a  final  term  involving  energy  absorption 
along  the  failure  boundary,  which  is  required  for  its  continued  growth.  It  is  pointed  out  that  this  latter 
term  can  result  in  severe  damping  of  the  radiation  field  and  is  therefore  of  importance  for  accurate 
representations  of  the  seismic  radiation  from  such  sources.  Wc  also  show  that  the  Green’s  integral  equa¬ 
tion  representing  the  seismic  radiation  field  can  be  expressed  in  equivalent  forms,  in  particular  in  a  form 
that  is  similar  to  that  obtained  earlier  using  an  initial  value  formulation  (eg.  Archambeau  and  Minster, 
1978)  and  also  in  the  form  of  a  combined  stress  pulse-dislocation  source  distributed  along  the  failure 
boundary.  One  advantage  of  the  present  formulation  is  that  it  precisely  prescribes  the  time  history  and 
spatial  variability  of  the  wave  fields  produced  in  terms  of  the  failure  zone  growth  rate,  geometry  and 
material  properties.  It  is  also  appropriate  for  the  general  case,  that  is  for  arbitrary  prestress  fields  and 
failure  growth  history  and  geometry.  Therefore,  a  moment  series  expansion  of  the  Green’s  integral 
representation,  based  on  the  directly  formulated  integral  equation,  or  any  of  its  equivalents,  will  produce 
moment  coefficients  that  are  analytically  related  to  the  fun  r  ;  tal  parameters  of  the  failure  process. 
Thus  the  results  of  this  study  provide  a  basis  for  inference  of  ' 'mice  physics,  as  well  as  a  physically 
based  predictive  capability  for  earthquake  sources. 
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I.  Introduction 


The  essence  of  a  kinematical  moment  tensor  representation  of  earthquake  seismic  radiation  (as 
opposed  to  a  dynamical  representation)  can  be  phrased  quite  simply.  In  particular,  the  source  is  con¬ 
sidered  to  be  representable  by  an  equivalent  body  force  term  in  the  equations  of  motion,  so  that  one 
assumes  that 

P  ^t2uk  ~  =  fk 

can  represent  the  source  effect  in  the  equation  of  motion,  with  fk  an  equivalent  body  force  distributed  over 
a  "source  volume"  and  of  a  form  which  will  give  rise  to  a  radiation  field  equivalent  to  that  of  the  actual 
source.  It  is  also  required,  in  effect,  that  the  "source  volume”  be  replaceable  by  an  elastic  region  in  which 
the  usual  linear  Green’s  function  representation  can  be  used.  In  this  context  a  dislocation  equivalent  to 
spontaneous  failure  in  a  stressed  solid  (that  is  failure  resulting  in  an  earthquake)  can  be  viewed  as  a 
choice  of  an  equivalent  body  force  that  is  distributed  over  a  planar  surface  within  the  medium.  (That  is, 
the  equivalent  body  force  has  a  delta  function  form  with  a  magnitude  that  is  usually  taken  to  be  propor¬ 
tional  to  an  imposed  displacement  offset  across  the  singular  "dislocation  plane”.)  Thus,  the  nonlinear 
effects  that  may  actually  occur  within  the  source  volume  are  simulated  by  an  equivalent  force  f  acting 
within  the  region  or  along  its  boundary.  In  this  case  the  dynamic  displacement  field  can  be  represented 
by: 

l 

um(r,t)  =  J  dt°  J  fk(r0,to)  Gkm(r,t;r0,to)  dr0 
n  o  v. 

where  Gk'(r,t;r0,t0)  represents  the  Green’s  tensor  describing  elastic  wave  propagation  in  the  medium  and 
V0  is  the  entire  volume  where  fk  is  non-zero  (See  for  example  Morse  and  Feshbach,  1953). 

An  example  of  an  expansion  for  u  in  terms  of  moments  of  f  may  be  obtained  by  formally  expanding 
the  Green’s  tensor  in  a  Taylor’s  series  in  the  source  coordinates  r0  (e.g.  Stump  and  Johnson,  1982),  so 


that: 


G“(r,t;r0,to)  =  ^ 


G  u(r^°-to) 


V<0)X(P).  .  .  _« 

\  xu  \ 


(0) 


where  summation  is  implied  for  the  repeated  coordinate  indices  1,, •  •  •  and  G$, ...  u  is  evaluated  at  r0  =  0 
(the  origin  of  the  source  coordinates),  with 


JU-U  1 


dxf  dxfj 


f Gk’(r,t;r0,t0)  | 

r,«0 


Using  this  expansion  for  G™  in  the  Green’s  function  representation  for  then  gives* 


um(r»t)  =  2<  ,  /  G)la1..!.(r,t;0,t0)  Mkli . . .  ^(to)  dto 
n‘  o 


where 

MkilJ.(t0)  =  J  xf  xf..xf>  fk(r0,to)  dr0 

v. 

is  the  (n  +  1)  th  order  moment  tensor. 

Generally  the  lowest  order  term  (n=  1),  corresponding  to  a  simple  point  double  couple  or  disloca¬ 
tion,  is  used  and  in  this  case: 

l 

ujr.t)  =  J  G£S(r,t;0,to)  Mu(to)  dto 

o 


Clearly,  all  of  this  is  completely  formal  with  no  actual  physical  description  of  the  source  involved. 
Further,  these  representations  follow  from  the  assumption  that  a  simple  distributed  equivalent  force,  act¬ 
ing  in  an  clastic  zone  replacing  the  true  physical  source,  can  actually  represent  (be  equivalent  to)  an  earth¬ 
quake  source.  That  this  is  not  necessarily  so  can  be  appreciated  by  reflecting  on  the  distinct  possibility 
that  the  true  equivalent  force  may  depend  explicitly  on  the  dynamic  displacement  field  itself  or  on  its 

+Thc  n  =  0  term  in  the  expansion  for  U  vanishes  by  the  conservation  of  momentum  when  there  are  no  external  forces  or 
torques  on  the  system. 
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space-time  derivatives.  In  addition,  in  the  Stump-Johnson  moment  expansion,  it  is  assumed  that  the 
series  expansion  for  G™,  and  resulting  moment  series  expansion  for  um,  exist  and  do  not  diverge.  While  it 
may  seem  intuitively  plausible  that  the  latter  assumptions  are  valid,  it  is  also  likely  that  expansions  of  this 
type  would  be  slowly  convergent,  and  certainly  extremely  so  for  a  complex  earthquake  of  large  dimen¬ 
sions. 

In  a  somewhat  similar  approach,  Archambeau  (1964,  1968)  considered  the  expansion  of  an 
equivalent  force  representation  in  moments,  except  that,  rather  than  formally  expanding  the  Green’s  func¬ 
tion  as  a  series,  the  equivalent  force  was  expanded  as  a  (vector)  harmonic  scries  and  the  Green’s  function 
integral  was  evaluated  as  a  moment  series  expansion.  While  this  approach  has  the  likelihood  of  produc¬ 
ing  a  more  rapidly  convergent  moment  series,  than  would  be  produced  by  the  expansion  of  the  Green’s 
function,  it  is  (as  was  pointed  out  in  the  original  development)  purely  formal  with  no  direct  relationship  to 
the  actual  source  physics. 

Therefore,  aside  from  these  mathematical  considerations,  the  essential  question  to  be  addressed  is: 
What  is  the  proper  equivalent  force  term  in  a  moment  expansion,  and  in  particular  is  there  an  equivalent 
force  that  can  be  related  to  the  actual  physical  source.  If  an  analytical  relationship  can  be  established, 
then  the  relationship  of  moment  tensor  observations  to  the  physics  of  the  source  can  be  established. 

As  will  be  shown,  the  complete  Green’s  function  representation  of  an  earthquake  source  results  in 
an  integral  equation  that  can  only  approximately  have  the  simple  equivalent  force  form.  But  even  assum¬ 
ing  the  earlier  formalism  to  be  approximately  valid  and  a  useful  representation  for  inversion,  before  the 
formal  inversion  results  can  properly  be  interpreted  it  is  still  necessary  to  be  able  to  relate  the  moment 
tensor  components  to  the  source  physics  in  a  rigorous  fashion. 

In  regard  to  the  question  of  the  physical  significance  of  the  moment  tensor.  Gilbert  (1970)  argued  on 
intuitive  grounds  that  the  equivalent  body  force  equivalent  could  be  expressed  in  terms  of  the  (static) 
stress  drop  resulting  from  an  earthquake  in  the  form 

fk  -  -  a,T;k 


3. 


where  Tfc  is  the  stress  change  (throughout  the  earth)  accompanying  the  seismic  event.  In  this  case,  the 
second  order  moment  tensor  is  just  T^.  Similarly,  Backus  and  Mulcahy  (1976a)  attempted  to  provide 
a  basis  for  the  choice  of  fk,  and  in  so  doing  introduced  the  idea  of  a  "stress  glut".  In  terms  of  the  stress 
glut  tensor,  this  equivalent  force  is 

fk  =  - 

where  r,*  denotes  the  "stress  glut".  For  an  earthquake  source,  r*  is  taken  as  the  difference  between  the 
true  physical  stress,  S^,  and  a  "model  stress".  with  this  latter  stress  function  corresponding,  within  the 
failure  zone,  to  the  spatial  continuation  of  the  linear  stress  from  outside  the  nonlinear  failure  zone  into 
this  zone.  Hence,  is  the  equivalent  of  a  "stress  drop"  within  the  failure  zone  and,  further,  the  second 
order  moment  tensor  is  equal  to  Tn.  in  this  case.  However,  r^,  by  virtue  of  its  definition,  is  zero  in  the 
surrounding  linear  zone,  so  that  while  Gilbert’s  equivalent  is  non-zero  everywhere,  the  stress  glut 
equivalent  is  non-zero  only  within  the  failure  zone  of  an  earthquake. 

In  this  paper  we  will  show  that  the  moment  tensor  is  actually  fundamentally  related  to  time  changes 
in  the  instantaneous  equilibrium  displacement  field  that  occurs  in  the  medium  surrounding  a  failure  zone, 
and  that  an  interpretation  of  the  moment  tensor  can  be  obtained  in  terms  of  changes  in  this  vector  field. 
Furthermore,  we  show  that  the  changes  in  the  equilibrium  displacements  can  be  related,  explicitly,  to 
failure  zone  shape,  growth  rate  and  material  rheological  properties,  as  well  as  to  the  initial  (tectonic) 
stress  in  the  medium.  This  relationship,  as  might  be  expected,  is  nonlinear  in  most  of  the  source  parame¬ 
ters,  so  that  inference  of  physical  parameters  from  the  moment  tensor  is  a  nonlinear  inversion  problem. 

In  order  to  develop  our  results  with  clarity  in  the  context  of  current  usage,  it  is  necessary  to  examine 
the  previous  moment  tensor  derivations  and  interpretations  in  some  detail.  Initially,  we  will  focus  our 
discussion  on  the  original  derivation  given  by  Gilbert  (1970),  and  the  critique  of  this  formulation  by 
Backus  and  Mulcahy  (1976a),  since  many  current  interpretations  derive  from  these  considerations.  In  this 
regard,  we  will  show  that  the  moment  tensor  cannot  be  related  to  a  "stress  drop"  in  the  manner  described 
by  Gilbert,  and  that  the  inclusion  of  gravitational  effects,  as  described  by  Backus  and  Mulcahy  in  their 
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critique  of  Gilbert’s  derivation,  docs  not  correct  the  problem.  In  fact,  the  changes  in  Gilbert’s  formula¬ 
tion  that  arc  advocated  by  Backus  and  Mulcahy  will  be  shown  to  lead  to  a  null  result;  that  is  the  source 
term  (and  therefore  the  seismic  radiation  field)  should  vanish  at  all  points  and  at  all  times  in  their  formal 
representation.  Gilbert’s  source  representation,  on  the  other  hand,  can  be  shown  to  only  include  the  very 
small  seismic  effects  due  to  gravity  and  density  changes  accompanying  an  earthquake. 

In  the  development  that  follows  in  the  next  section  we  will  first  show  why  Gilbert’s  phcnomcnolog 
ical  moment  tensor  formulation  is  incorrect  and  why  the  argument  involving  inclusion  of  gravitational 
effects,  as  advanced  by  Backus  and  Mulcahy,  is  similarly  incorrect.  We  will  also  show  that  Backus  and 
Mulcahy’s  moment  tensor  representation,  including  the  ‘stress  glut’  representation,  docs  not  properly 
represent  the  physics  of  a  spontaneous  seismic  source,  such  as  an  earthquake,  nor  docs  it  properly 
represent  an  earthquake  radiation  field  in  an  equivalent  sense.  We  will  then  demonstrate  that  a  basically 
different  approach  is  required  in  order  to  properly  describe  the  physics  of  a  spontaneous  seismic  source. 
In  particular,  we  will  show  how  a  spontaneous  failure  process  is  properly  formulated  dynamically  and, 
when  so  described,  provides  a  non-trivial  moment  tensor  that  can  be  related  to  the  physics  of  the  source. 

II.  Gilbert’s  "Phenomenological"  Moment  Tensor  Representation 

Gilbert’s  (1970)  original  development  of  a  moment  tensor  representation  followed  the  essentially 
intuitive  approach  that  is  still  commonly  evoked  to  model  an  earthquake;  namely,  it  was  assumed* 

(i)  That  the  external  body  force  term  f  in  the  equations  of  motion: 

P  d,2uk  -  <5itu  =  p  fk 

can  be  used  to  represent  the  (force)  effects  associated  with  an  earthquake. 

(ii)  That  the  initial  state  of  the  medium  (at  t  =  0)  is  described  by: 

uk(r,0)  =  3luk(r,0)  =  0 

* - 

Gilbert  used  equations  of  motion  for  a  discrete  system  of  panicles  and  passed  to  the  limit  of  a  continuum  for  his  final 
results  In  this  discussion  continuum  equations  are  used  from  the  onset  since  use  of  the  discrete  system  is  neither  neces¬ 
sary  nor  particularly  useful. 


(iii)  That  the  earthquake  may  be  adequately  represented,  for  the  purposes  of  predicting  the  normal 
mode  response  of  the  planet,  by  a  step  function  response;  so  in  view  of  the  first  assumption: 

fk(r.l)  =  Fk(r)  H(l) 


(iv)  That  ‘the  body  force  caused  by  the  stress  drop  is’: 

Fk(r)  =  -  3iTu 

where  the  spatial  stress  drop  is  defined  in  terms  of  the  stress  after  an  earthquake,  xfip,  and  the 
stress  before  an  earthquake,  x£\  as: 

Ttt(r)=  f-x^jtX) 

These  assumptions  appear  eminently  plausible  at  first  sight.  However,  they  lead  directly  to  a  con¬ 
tradiction  which  arises  from  the  fact  that  both  x£J  and  x^  satisfy  the  equations  of  equilibrium.  That  is: 


=  -  P  gk®  :  r  £  V0  . 

0) 

fltuX J  =  0 ;  r  e  3V0 

where  the  bracket  notation  applied  to  x^n,  on  3V0,  the  boundary  of  V0,  denotes  the  change  in  the  traction 
across  the  boundary.  Likewise: 


3|tf)*-pgf);reVl 
=  0 ;  r  e  3V, 


(2) 


where  pgkm  and  pgk(F)  represent  the  gravitational  forces  acting  and  where  the  boundary  conditions  express 

continuity  of  tractions  at  (all)  the  medium  boundaries.  Here  the  field  x,^  is  defined  over  the  volume  V0 

consisting  of  the  entire  volume  of  the  planet,  with  3V0  representing  its  boundaries,  while  x^p  is  defined 
_ _  “ 

In  addition  to  the  free  surface  of  the  planet,  there  arc  other  boundaries  within  the  earth  where  the  elastic  properties  may 
change  disconlinuously  and  which,  therefore,  must  be  accounted  for  as  internal  boundaries  in  Vq  and  Vj.  However,  these 
can  be  omitted  from  explicit  mention  if  we  observe  once  and  for  all  that  the  tractions  are  continuous  at  all  such  boundaries 
and  that  at  solid-solid  boundaries  the  displacement  is  continuous,  while  on  fluid-solid  or  fluid-fluid  boundaries  only  the 
normal  component  of  the  velocity  is  continuous.  In  the  present  context  these  boundaries  and  the  conditions  on  them  do  not 
change  (to  first  order)  before  and  after  an  earthquake,  so  the  same  conditions  will  always  apply  and  are  to  be  assumed 
throughout. 
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over  V[  consisting  of  the  ntcrior  of  the  planet  outside  the  failure  /one  of  the  earthquake,  and  having 
boundaries  denoted  by  3V,,  which  includes  the  (newly  formed)  failure  /one  boundary  in  addition  to  the 
boundaries  3V0,  existing  prior  to  the  event.’* 

Inside  the  failure  zone  we  have 

d,T,k  =  -p'fv ;  re  V'  =  V0OV, 
hihk  ~  hk  >nk  ;r  edv' 

where  i,k  is  any  (residual)  stress  which  may  exist  withn  the  failure  zone  after  the  earthquake,  while  is 
the  final  equilibrium  stress  in  the  clastic  zone  surrounding  the  failure  zone  and  3V’  denotes  the  failure 
zone  boundary.  Here  the  symbol  ©denotes  the  set  theoretic  difference.  The  regions  V0,  V,  and  V'  and 
their  boundaries  arc  schematically  illustrated  in  Figure  1. 

It  is  important  to  emphasize  the  necessity  of  taking  explicit  account  of  the  existence  of  the  failure 
zone  and  to  express  the  equilibrium  equations  in  each  zone  separately,  with  the  appropriate  boundary  con¬ 
ditions  serving  to  connect  the  stress  fields  within  the  non-linear  source  region  to  those  in  the  linear  region 
surrounding  it  (i.c.  the  rest  of  the  earth).  This  necessity  arises  from  two  tacts.  First,  the  only  change  in 
the  boundaries  within  the  earth  is  the  formation  of  a  failure  zone  boundary  and,  while  there  may  be  many 
other  boundaries  along  which  material  discontinuities  occur  and  where  boundary  conditions  apply,  this  is 
the  only  boundary  at  which  non-negligible  changes  occur  during  an  earthquake.  Second,  and  most  criti¬ 
cal  for  the  ensuing  discussion,  the  failure  zone  defines  a  region  of  non-linear  behavior  during  its  forma¬ 
tion,  within  which  it  is  not  possible  to  apply  the  usual  linear  theory.  As  will  be  seen,  part  of  the  problem 
with  Gilbert’s  representation  arises  from  the  fact  that  he  does  not  account  for  the  existence  of  the  failure 
zone  boundary  and  the  associated  boundary  conditions  and,  in  effect,  assumes  linearity  everywhere 
throughout  the  planet  at  all  times,  including  within  the  failure  /one.  This  particular  problem  also  arises, 
in  .  more  direct  fashion,  in  the  equivalent  source  representation  advocated  by  Backus  and  Mulcahy 
(1976a, b)  and  will  be  discussed  later. 
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In  order  to  relate  the  description  of  equilibrium  fields  given  in  (1)  through  (3)  to  Gilbert’s  use  of 
equilibrium  fields  defined  on  the  interior  of  the  planet,  we  observe  that  the  region  over  which  Gilbert 
defines  the  stress  drop,  in  (iv),  is  denoted  as  V  by  him  and  corresponds  to  our  region  V0  (i.e.  the  entire 
planet).  To  properly  define  a  stress  drop  over  the  entire  interior  of  the  planet  it  is  necessary  to  define  it 
as: 

Tlk  =  t^)-t1f  ;reV, 

r  e  V  =  V0©V 

and  to  include  the  boundary  conditions  applicable  on  the  failure  zone  boundary  dv ',  in  order  to  connect 
the  stress  drop  in  V  ’  with  that  outside  the  failure  zone,  in  V,.  This  condition  is  simply, 

Tu,n,  =  ;  r  e  3V '  (4-b) 

which  equates  the  tractions  at  the  boundary.  Thus,  while  the  tractions  are  continuous,  the  stress  drop 
components  themselves  may  be  discontinuous  across  3V '. 

Now  we  can  relate  this  stress  drop  definition  to  Gilbert’s  use  of  stress  drop  by  simply  setting  the 
defined  in  (iv)  to  be  equal  to  one  or  the  other  of  the  expressions  in  (4-a),  depending  on  whether  the  coor¬ 
dinate  point  is  in  V,  or  V This  modification  of  Gilbert’s  definition  can  be  thought  of  as  an  extension  of 
his  definition,  required  in  order  that  the  failure  zone  region  is  explicitly  (and  correctly)  "covered"  with  the 
correct  boundary  conditions  included.  We  shall  show  below,  however,  that  this  is  not  the  main  difficulty 
with  Gilbert’s  forms!  rcr  ills,  even  though  a  proper  definition  of  the  failure  zone  equilibrium  and  boun¬ 
dary  condition1  >  jndamcntal  to  the  correct  formulation  of  the  problem. 

In  this  regard,  vj  note  from  (1)  and  (2)  that  while  and  x IP  are  themselves  different,  because 
they  satisfy  different  set  of  boundary  conditions  (one  with  the  failure  zone  present,  the  other  without  it), 
the  divergences  of  these  stresses  can  be  related,  since  from  (1)  and  (2); 

dijxiiP  -  Tifj  =  -  •  p  gk0)  -  P  gk(F)j  J  r  £  v, 
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by  simple  subtraction  of  the  equations  of  equilibrium.  As  indicated,  the  relation  applies  over  the  region 
exterior  to  the  failure  zone  (V  ).  Now  we  sec  from  Gilbert’s  assumption  (iv)  that  the  quantity  in  the 
brackets  is  defined  as  the  stress  drop,  and  so  the  stress  drop  must  satisfy  the  relation: 


Here  pg^  represents  the  new  gravitational  force  field  acting  after  the  earthquake,  whereas  pgk0)  is  the  field 
before  the  event.  The  difference  in  these  force  fields  is  due  to  llic  redistribution  of  mass  due  to  the  earth¬ 
quake,  resulting  in  changes  in  density  and  associated  changes  in  the  gravitational  acceleration.  Gilbert 
docs  not  directly  concern  himself  with  gravitational  effects.  (He  docs  not  neglect  gravity  effects  as 
claimed  by  Backus  and  Mulcahy  (1976a),  he  simply  asserts  that  the  body  force  equivalent  of  an  earth¬ 
quake  is  -3,1^  and  never  observes  that  (5)  is  actually  true.)  However,  if  gravity  effects  are,  in  fact, 
neglected  in  comparison  with  effects  of  tectonic  origin,  as  is  shown  to  be  justified  in  the  Appendix  1, 
then: 

a1Tlk  =  0;reV1 .  (6-a) 

We  also  observe  that  the  stress  drop  within  the  non-linear  failure  zone  is  also  equal  to  the  gravity 
field  change,  that  is: 

=  -  Pgkn)-P8k]  :  reV', 

and  when  gravity  changes  arc  ignored  relative  to  tectonic  effects,  then  we  also  have 

a,T»  =  0  ;  r  e  V  '  (6-b) 

within  the  failure  zone. 

Therefore  we  sec  from  (6)  that,  to  the  approximation  in  which  neglect  of  the  gravitational  forces  is 
appropriate,  the  divergence  of  the  stress  drop  is  zero  everywhere.  In  terms  of  Gilbert’s  (equivalent) 
‘body  force  caused  by  the  earthquake’,  as  defined  in  (iv).  we  have  that  this  quantity  is  also  zero.  That  is, 
we  have: 
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Fk(r)s-9!T,k  =  0;reV0, 


and,  ae  well, 

fk(r,t)  =  Fk(r)  H(l)  =  0 ; r  e  V0 , 

in  view  of  (6). 

Therefore  the  body  force  term  used  by  Gilbert  to  derive  the  moment  tensor  representation  is,  in  fact, 
zero  (or  at  least  negligible)  and  the  result  derived  from  the  normal  mode  representation  using  this 
equivalent  force  does  not  represent  an  earthquake  displacement  field.  In  particular,  Gilbert  finds  that  the 
displacement  field  due  to  a  body  force  term  of  the  form  fk(r,t)  =  F,(r)  H(t),  with  H(t)  a  step  function,  can  be 
put  in  the  form 


where  the  s(n)  arc  the  normal  mode  displacement  eigenfunctions  for  the  entire  "unperturbed"  planet;  that  is 
the  planet  without  the  failure  zone  boundary  present  within  it.  Here  s<n)  are  the  complex  conjugates  of  the 
mode  functions,  are  the  mode  angular  frequencies  and  Q„  the  associated  dissipation  function  for  each 
mode.  (For  this  result,  the  Qn  are  assumed  to  be  much  larger  than  unity  so  that  dissipation  is  small.)  This 
expansion  is  formally  correct  if  effects  from  the  failure  zone  boundary  arc  neglected  and  if  it  is  applied  in 
the  region  V,  (outside  the  nonlinear  failure  zone). 

However,  insertion  of  -  3^  for  Fk  in  this  representation,  whether  intuitively  plausible  or  not,  is  ad 
hoc  and  quite  incorrect.  In  fact,  as  we  have  noted,  if  gravity  effects  arc  negligible  relative  to  tectonic 
effects  (which  is  in  fact  true),  then  the  divergence  of  the  stress  drop  vanishes  and  so  does  the  right  side  of 
(7).  Therefore,  Gilbert’s  application  of  Gauss’  theorem  to  the  volume  integral  term  involving  the  diver¬ 
gence  of  the  stress  drop  and  subsequent  approximations  which  lead  to  the  widely  used  result  (e.g.  Gcllcr, 
1976;  McCowan  and  D/.icwonski,  1977;  Dzicwonski  and  Gilbert,  1974;  Gilbert,  1973;  Gilbert  and 
D/.icwonski,  1975;  McCowan,  1976;  Luh  and  Dzicwonski,  1976) ; 
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{m*  e<" j 

1  -  cosoy  exp 

k‘j 

Gin 

. 

(X) 


is  without  proper  foundation  and  the  result  is  incorrect.  (Here  Mu  is  the  "moment  tensor",  which  in  this 
result  equals  the  stress  T^,  while  e^n)  is  the  strain  tensor  associated  with  the  complex  conjugate  displace¬ 
ment  eigenfunctions  for  the  planet.)  Basically,  the  error  arises  from  the  fact  that  the  divergence  of  the 
stress  drop  is  not  a  proper  equivalent  body  force  for  an  earthquake  source,  plausible  ad  hoc  assertions  to 
the  contrary.  The  fact  that  this  equivalent  force  actually  vanishes  when  gravity  effects  arc  neglected  is 
simply  indicative  of  the  impropriety  of  the  assumption. 


The  expression  in  (8)  is  not,  however,  zero  as  would  be  expected.  That  is,  since  Fk  =  -  3,^  vanishes 
(to  the  approximation  in  which  gravity  changes  are  neglected)  and  so  causes  (7)  to  give  a  null  displace¬ 
ment  field  result,  one  might  certainly  expect  that  any  result  derived  from  (7),  such  as  (8),  should  also 
vanish.  This  is  not  the  case  because,  in  deriving  (8)  from  equation  (7),  Gilbert  omits  the  surface  integral 
term  over  the  failure  zone  boundary  so  that  (8)  represents  only  part  of  the  solution  (7).  The  part  neglected 
(the  surface  integral  term)  exactly  cancels  the  right  side  of  (8),  to  give  the  (correct)  null  result  or,  if  grav¬ 
ity  is  not  neglected,  the  displacement  effects  due  to  gravity  field  changes.  That  is,  if  we  directly  substi¬ 
tute  the  divergence  of  the  stress  drop  into  Gilbert’s  equation  (7),  then  we  get: 


um(r,t)  =  Xsin)(r>  J 


fP  E™  ~  P  g^l  Skn)dV0 


1  -  cos  iont  cxp(-  (ont/2Qn) 


(9) 


since  the  divergence  of  the  stress  drop  has  the  form 

diTtt  =  [p  gk(I)  -  P  gkao 

everywhere  :n  the  planet,  including  inside  the  failure  zone.  Thus,  when  the  divergence  of  the  static  stress 
drop  is  used  to  represent  the  "equivalent  force"  due  a  spontaneous  source  like  an  earthquake,  then  the  dis¬ 
placement  predicted  will  be  negligibly  small  or  zero  and  not  representative  of  the  seismic  radiation  field 
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produced  by  the  event. 


Backus  and  Mulcahy  (1976a)  on  the  other  hand  argue  that  Gilbert’s  results  are  in  error  precisely 
because  of  the  neglect  of  gravity  in  defining  a  moment  tensor.  But  this  is  not  the  reason  that  the  result  is 
in  error,  as  was  already  indicated  above.  In  particular,  Gilbert  does  not  neglect  gravity  since  he  simply 
assumes  that  the  divergence  of  the  stress  drop  is  the  equivalent  force  to  be  associated  with  an  earthquake 
and,  once  this  is  done,  gravity  effects  arc  actually  implicitly  included.  In  fact,  as  we  have  shown  in  equa¬ 
tion  (9),  the  true  consequence  of  Gilbert’s  assumption  is  that  he  actually  should  obtain  only  that  part  of 
the  displacement  due  to  an  earthquake  that  is  associated  with  gravitational  field  changes.  Further,  as  will 
be  directly  demonstrated,  the  formulation  suggested  by  Backus  and  Mulcahy  gives  a  null  result  when 
gravity  is  included  in  the  way  they  advocate. 

In  particular,  if  the  gravitational  changes  arc  included,  as  suggested  by  Backus  and  Mulcahy  (1976a, 
p.  346),  by  defining  a  gravitational  stress  tensor  (Gy)  such  that 

P  Sj  ^iGy  , 

where 

Gy  =  (87tG)"’  (gkgk8;j  -  2gigJ) , 

then  the  equilibrium  equations  for  the  initial  and  final  states,  equations  (1),  (2)  and  (3),  have  the  form: 

difx^P  +  G£>]  »  0 ;  r  e  V0 

d.[^;,  +  Or]=0;re  V,  .  .  (10) 

di  %  +  Gi]  =  0 ;  r  E  V0  0  V, 

The  formal  expressions  given  by  Backus  and  Mulcahy  (1976a,  p.  347)  for  the  equilibrium  equations  arc, 
however  (in  our  notaiion):+ 

+Uackus  and  Mulcahy  use  the  term  "all  space",  which  is  apparently  meant  to  be  taken  literally.  We  shall  confine  our  dis¬ 
cussion  somewhat  more  to  the  point  however,  and  take  the  liberty  of  replacing  their  "all  space"  by  the  entire  interior  of  the 
earth.  In  our  notation  this  is  V0. 
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3i  T]jP  +  G^j  -  0  ;  r  e  Vo 


3]  if  +  GiiPj  =  -  Yk  :  r  e  V0 

where  represents  ‘the  equivalent  volume  distribution  of  body  forces’  and  is  asserted  to  be  Gilbert’s 
‘equivalent  body  force  due  to  the  source,  corrected  for  gravitational  effects’. 

Here,  since  Backus  and  Mulcahy  ignore  the  boundary  conditions  on  the  failure  zone,  as  did  Gilbert, 
and  define  the  equilibrium  equation  after  the  failure  process  over  the  entire  interior  of  the  planet,  then  we 
must  define  the  stress  tensor  for  the  final  equilibrium  state  to  be  Tp  when  r  e  V,  (inside  the  clastic  zone  of 
the  planet)  and  to  be  %  when  re  Vu©V,  (inside  the  failure  zone),  as  was  done  earlier  in  (4-a).  Thus,  if 
we  write  Backus  and  Mulcahy’s  expressions  over  the  range  of  r  for  which  they  define  it,  but  separately 
for  r  in  V,  and  .c :  r  in  V  (the  failure  zone),  we  have: 

3,  tP  +  GP]  =  re  V, 

3i[%  +  Gi]  re  V' 

However,  comparing  these  directly  with  the  required  equilibrium  equations  given  in  (10)  shows  that: 

Yk  =  0 ;  r  e  V , 

y{  =  Q;reV' . 

That  is,  as  must  be  true,  equilibrium  requires  to  vanish  everywhere,  including  within  the  failure  zone 
V  Therefore  it  is  clear  that: 

Ykv  =  0;reV0=  V,©V' 

so  that  this  body  force  equivalent  is,  in  fact,  identically  zero  in  both  V,  and  V  =V0©V,  (that  is, 
throughout  the  planet  including  the  region  within  the  failure  zone).  Thus,  it  contributes  precisely  nothing 
to  the  radiation  field  from  an  earthquake  for  essentially  the  same  reason  that  Gilbert's  stress  drop 
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equivalent  force  gives  a  null  result  when  gravity  is  neglected  in  his  formulation.  Thus,  use  of  the  diver¬ 
gence  of  changes  in  the  equilibrium  stress-gravity  field  as  the  driving  force  for  a  spontaneous  source  will 
necessarily  give  a  null  result,  since  the  mathematical  expression  for  the  existence  of  equilibrium  is  pre¬ 
cisely  that  this  particular  force  vanish. 

As  a  consequence  of  these  results,  Backus  and  Mulcahy’s  definition  of  an  extended  moment  tensor 
density  as 

Mu  =  (T,Sr)-'t£))  +  (Gff)-Gk<?) 

and  subsequent  arguments  about  the  size  of  the  gravitational  term  (G^P  -  G^)  relative  to  the  stress  drop 
term  (xjP  -  t^)  are  moot,  since 

diMu  =  -  Yk  ■  0 

in  any  case.  Thus,  since  the  divergence  of  Mu  is  the  actual  source  term,  we  see  that  this  "source  term" 
must  necessarily  vanish,  whatever  the  magnitudes  of  the  stress  drop  and  gravity  field  changes. 

Backus  and  Mulcahy  also  define  the  boundary  conditions  applicable  to  the  equilibrium  states  before 
and  after  the  earthquake,  and  state  them  as  (in  our  notation): 

n,  t,$P  +  G|$P  =  0 ;  on  3  V 
n,  +  on  3V 

where  3V  denotes  the  outer  boundary  of  the  planet  (surface  of  the  earth),  since  V  is  used  by  them  to 
denote  the  whole  interior  region  of  the  earth.  Here  jump  notation  is  used,  so  the  plus  and  minus  signs  on 
the  brackets  denote  the  difference  in  the  bracketed  quantity  across  the  boundary  3V,  when  approached  in 
the  limit  from  opposite  sides.  As  noted  earlier  they  ignore  the  failure  zone  boundary  condition,  since  dV 
(which  is  equivalent  to  our  boundary  3V0)  is  the  surface  of  the  earth.  They  introduce  the  surface  source 
term  yk’  to  include  the  case  where  the  failure  zone  intersects  the  surface  of  the  earth,  in  which  case  it  is 
claimed  that  yk  will  give  the  proper  equivalent  force  term  along  this  intersection  on  3V. 
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However,  inclusion  of  this  "surface  distributed"  source  term,  contributes  nothing  of  substance  to 
their  formulation,  and  it  still  gives  a  null  result.  That  is,  consider  the  ease  when  the  failure  zone  docs  not 
happen  to  intersect  the  earth’s  surface.  Then  =  0  by  definition.  Clearly  then  all  possible  equivalent 
source  terms  vanish,  that  is  =  0  and  ^0,  and  Backus  and  Mulcahy’s  formulation  gives  a  null  result. 
Thus,  it  predicts  no  displacement  field  resulting  from  the  event  whatsoever,  because  all  source  terms  van¬ 
ish,  and  this  is  obviously  incorrect.  Since  this  reformulation  of  Gilbert’s  representation  is  incorrect  for  all 
earthquakes  for  which  the  failure  zone  docs  not  intersect  the  earth’s  surface,  it  is  reasonable  to  conclude, 
without  going  through  a  detailed  evaluation  of  the  special  ease  of  a  failure  zone  intersecting  the  free  sur¬ 
face  of  the  earth,  that  the  reformulation  is  also  incorrect  in  general.  (In  fact,  since  Yk  -  0  everywhere  in  V, 
then  it  can  be  shown  that  Yk  must  also  be  zero  on  3V.) 

At  this  point  we  have  rather  thoroughly  explored  die  so-called  stress  drop  moment  tensor  represen¬ 
tation  of  Gilbert  and  the  ‘extended  moment  tensor’  representation  of  Backus  and  Mulcahy,  which  is 
intended  to  correct  Gilbert’s  formulation.  However,  we  have  shown  both  to  be  misdirected  and  to  give 
incorrect  representations  of  an  earthquake  source  radiation  field.  Indeed  both  give  (essentially)  null 
results  when  the  required  equilibrium  equations  are  taken  into  account.  Thus,  neither  of  these  formula¬ 
tions  account  for  the  physics  of  a  spontaneous  failure  source,  nor  do  they  kinematically  or  phenomenolog¬ 
ically  represent  such  a  source. 

III.  The  Backus-Muicahy  Stress  Glut  Formulation 

The  main  point  of  the  work  by  Backus  and  Mulcahy  (1976a,b)  is  not  the  reformulation  of  Gilbert’s 
stress  drop  equivalent,  but  the  formulation  of  a  comprehensive  phenomenological  moment  tensor 
representation  for  seismic  sources  in  general  and  earthquake  sources  in  particular.  It  is  therefore 
appropriate  to  comment  upon  the  basis  of  their  representation  since  we  intend  to  generate  a  type  of 
moment  tensor  representation  in  a  following  section,  but  from  quite  a  different  point  of  view.  Further, 
their  results  arc  very  similar  to  Gilbert’s  and  arc  similarly  incorrect. 
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The  essence  of  Backus  and  Mulcahy’s  approach  (Backus  and  Mulcahy,  1976a;  p.  342-346)  is  based 


on  the  following  considerations  (in  their  notation): 

(i)  The  exact  equations  of  motion  and  gravity  for  the  earth  are 

P  a;  =  a.S.j  -  p  +  fjv  in  V 

n,Sy  =  f*  on  3  V 
V2^  =  4nGp  ,  everywhere 
V  and  3i\jr  continuous  on  3V 
-» 0 ,  at  infinity 

where  V  is  the  (entire)  volume  of  the  earth,  3V  its  surface  and  fv,  f*  are  externally  applied  body 
and  surface  forces  respectively.  Here  y  is  the  gravitational  potential  at  any  point,  so  gj  =  3^  is 
the  gravitational  acceleration.  Finally,  S,,  are  the  components  of  the  'true  physical  stress'  in  the 
earth  and  aj  is  the  particle  acceleration. 

(ii)  A  'Mathematical  model  stress’  may  be  defined  as: 

Cy =  Eijiu  3ks, 

where  Eijkl  is  the  ordinary  elastic  tensor,  or  as: 

Cj =  Sj  -  sk3kS°  +  Fijk]  3kSj 

when  prestress  S°  and  gravity  arc  taken  into  account/  Here  sk  is  the  (seismic)  displacement  in 
the  elastic  region  of  the  earth.  In  addition,  Fjjld  is  the  elastic  tensor  corrected  for  the  presence  of 
prestress,  so  that.-++ 

+The  term  Sg3gS  ®  is  always  small  relative  to  the  other  terms  in  in  scismological  applications,  and  can  always  be 
neglected. 

^Since  real  stresses  in  the  earth  are  orders  of  magnitude  less  than  the  magnitudes  of  the  elastic  constants  appearing  in 
Eijkl,  then  it  is  justified  to  take  Fjjgj  =  E,^  in  most  scismological  applications  and,  in  particular,  in  this  one. 
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F,jki  =  E^i  +  1/2 -  Sv.S"  +  5,jS^  -  5;kS,;  +  5^  -  S^S" 

(iii)  The  model  stress  can  be  calculated  from  the  true  displacement,  sk,  by  virtue  of  its  definition  and, 
with  the  introduction  of  the  ‘model  stress’  in  the  equations  of  motion  (which  can  be  achieved 
by  simply  adding  the  divergence  of  this  ‘model’  stress  field  to  both  sides  of  the  equation  of 
motion  in  (i)),  one  obtains  equations  of  motion  expressed  in  terms  of  the  model  stress  plus  an 
extra  ‘source  term’,  which  is  used  by  Backus  and  Mulcahy  as  the  "equivalent  driving  function" 
or  "equivalent  source  function"  for  a  spontaneous  release  of  energy  corresponding  to  an  earth¬ 
quake.  In  particular,  the  equations  of  motion  become,  from  (i): 

P  +  d.Cij  =  d.Ci,  +  3,S0  -  P  fyl/  +  f/  .  in  V, 

when  the  divergence  of  the  model  stress  is  (formally)  added  to  both  sides  of  the  equation  of 
motion  in  (i).  This  equation  can  also  be  written  as: 

paj=aiciJ-pajv+f;+Y;,inv, 

with, 

tf  =  a.sa-a& 

where  y)v  is  viewed  (by  Backus  and  Mulcahy)  as  an  equivalent  source  term  arising  from  the 
introduction  of  the  model  stress  Cij-  Further,  the  boundary  condition  on  3V  then  becomes,  in 
terms  of  the  model  stress: 

niCij -  f/  +  Yj‘  .on  3V 

where 

7?  =  n,Cij  -  n.Sjj ,  on  dV 

with  Yj8  considered  to  be  (by  Backus  and  Mulcahy)  a  boundary  source  term,  on  av,  at  the  surface 
of  the  planet. 
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(iv)  A  moment  tensor  density  M(J  is  defined  in  V  such  that 

in  V, 

n,  =  7)* ,  on  9V. 

In  addition,  the  "stress  glut"  r.j  is  defined  to  be  the  difference  between  the  model  stress  ^  and 
the  true  physical  stress  S,j,  so: 

Hj  =  Cij  _  Sy  . 

Backus  and  Mulcahy  assert  that  'evidently  the  physical  source  region  is  precisely  the  region 
where  the  stress  glut  is  non-zero’  This  statement  therefore  fixes  the  choice  for  the  model  stress 
to  be  the  (essentially)  elastic  stress  outside  the  failure  zone,  which  is  identical  to  their  ‘true 
physical  stress’  in  the  elastic  zone. 

(v)  The  equations  of  motion  and  gravity,  given  in  (i)  above,  are  linearized  throughout  the  volume 
V  (the  entire  planet  including  the  failure  zone  volume)  and  a  relative  model  stress  £,■  is  intro¬ 
duced  such  that: 

where  S°  is  the  prcstrcss  in  V.  The  equations  of  motion  and  gravity  are  then  expressed  as+: 
p°  ay  =  a&j  -  p°  ay  -  P*  ay  +  r;  +  ,  in  v, 

n,  C,j  +  ak(s kS°)  -  (aks,)  Sk°]  =  fj*  +  Yj‘ ,  on  av, 

p1  -  -  ak(p°  sk) ,  inside  and  outside  of  V, 


VV  =  4nGp'  ,  inside  and  outside  of  V, 

_  « 

+As  noted  earlier  in 
the  prestress  Sy  and 

non-linear  zone,  but  these  linearized  equations  are  not  valid  in  this  zone  in  any  case. 


(ii),  the  various  terms  in  the  boundary  conditions,  stress  tensor  and  elastic  tensor  Fljkj,  involving 
derivatives  of  it,  are  small  in  the  clastic  zone  and  ordinarily  are  neglected.  TTiey  could  be  large  in  the 
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\y'j_  =0  and  n,  3k\)/!  =4nGplln,s1  ,  on  r)V 


\j/1  — >  0  ,  at  infinity, 


Cj  =  -  SiAS.j  +  Fljklf)kSi ,  inV. 

In  addition,  the  initial  conditions 

Sj  =  9,-Sj  =  0 ,  in  V  at  t  -  0, 


are  used. 

A  formal  eigenfunction  expansion,  using  the  norma!  modes  of  the  earth,  is  then  used  to  express  a 
solution  of  these  linear  equations  as: 


s/x.t)  =  £  Av(t)  u/v)(x) 

V=1 


with 


J  y/  Uj(v)  dv  +  J  y,5  Uj(v)  da  4  J  r;  o,M  d  v  + J  r;  u/v)  d<al 

v  av  v  av 


where  the  integration  is  over  the  entire  volume  V  and  outer  surface  dV  of  the  planet  and  u/v>  arc  the  spa¬ 
tial  eigenfunctions  for  the  planet,  with  denoting  the  complex  conjugate/ 

The  Backus  and  Mulcahy  result  for  an  earthquake  source,  where  no  external  body  forces  arc  applied 
and  the  ’source  zone’  docs  not  intersect  the  earth’s  surface  dv,  is  simply 


dt  sin  o\,(t  -t)  J  (a.r.j)  U/v)  dV 

v 


*Both  Gilbert  (1970)  and  Backus  and  Mulcahy  (1976a)  cite  completeness  of  the  modal  eigenfunction  set  as  the  basis  for 
their  representations  of  the  seismic  field  from  an  earthquake.  But  since  they  completely  omit  the  boundary  condition  at  the 
failure  none,  they  end  up  omitting  integral  terms  in  Av(t),  the  "excitation  function",  that  represent  scattering  from  the 
failure  zone  boundary,  among  other  effects.  (Such  integral  terms  arc  explicitly  considered  in  lalcr  sections  of  this  paper  and 
are  shown  to  ansc  quite  naturally  when  a  proper  Green's  function  integral  representation  is  used  to  describe  the  seismic  ra¬ 
diation  effects  associated  with  a  spontaneous  failure  process  -  that  is  an  earthquake.)  Further,  and  most  important,  these  au- 


ukv)(x) 


01) 


Av(t) 


I 

=  w;1  J  dt 


sinco^l-  t)i 


19 


in  terms  of  the  ‘stress-glut’.  Here  the  spatial  integration  is  only  over  the  ‘source  region’,  that  is  V  the 
failure  zone,  and  the  ‘model  stress’  is  the  stress  in  the  linear  region.  (If  the  ‘source  region’  intersects  the 
free  surface  3V,  then  their  result  would  also  include  the  surface  integral  over  3V  involving  Y|*  =  n.r, ,.) 

The  logical  development  described  in  (i)  through  (v)  is  not,  however,  a  proper  one  and  the  ‘stress 
glut’  and  moment  tensor  representations  involving  volume  integration  over  the  failure  zone  are,  conse¬ 
quently,  incorrect. 

The  precise  nature  of  the  problem  with  their  approach  can  be  seen  most  clearly  by  considering  the 
failure  zone,  and  its  associated  boundary,  explicitly  in  the  formulation.  We  employ  the  notation  used  ear¬ 
lier  and  indicated  in  Figure  1.  Appropriate  equations  of  motion  in  the  separate  regions  V  (the  nonlinear 
failure  zone)  and  V,  (the  linear  zone  surrounding  the  failure  zone)  are 

^  [pvt]  ~ djS.j  =  p'g,'  ;  r  e  V 
PAV,0)  -  3jS,f  =  p,g,(1)  ;  reV, 

with  the  tractions,  S.'n,  and  S^nj,  ind  velocity  fields,  v  and  v(1),  in  V '  and  Vj  being  related  by  boundary 
equations  expressing  conservation  of  momentum,  mass  and  energy  across  dv ',  the  failure  zone  boundary. 
On  the  other  hand  we  see  from  (i)  above  that  Backus  and  Mulcahy  write  the  equations  of  motion 
throughout  the  planet  in  lire  form 

P  a,  -  =  pg,  (13) 

with  a,  the  acceleration.  (Applied  external  forces  arc  not  relevant  and  arc  omitted  throughout  the  remain¬ 
ing  discussion.)  This  equation  can  be  formally  related  to  those  in  V  and  V,  by  simply  letting  the  equation 
take  the  form  of  one  or  the  other  of  the  equations  in  (12),  depending  on  whether  r  is  in  V  or  V,.  Since 
Backus  and  Mulcahy  use  the  stress  in  the  linear  zone,  S^,  as  their  model  stress  (which,  in  fact,  is  the  only 
choice)  then,  as  in  (iii)  above,  this  gives  for  their  equation  (13): 

thors  consequently  ignore  the  essential  nonlinear  properties  of  the  failure  /an c  and  integrate  through  this  region,  thereby 
summing  (physical  or  "lme"j  stress  changes  and  gravity  field  changes  from  within  this  nonlinear  zone  as  if  they  were  linear 
effects  which  superpose  linearly  and  are  propagated  linearly  within  this  zone. 
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pa1-3JS1<J,)  =  pg,  f  y,v  (14) 

with 

Y|V  =  3j(Sy  ~  S,j')) 

Now  in  terms  of  the  equations  given  in  (12)  for  the  linear  and  non-linear  regions,  this  referencing  to  S,'1’ 
simply  gives: 


dt 


Pid<vl(l)-ajSiy,  =  p1g1U)  ;  re  V, 


p'vij  -  djS^  =  pg,'  +  3j[S,j  -  Sf|  ;  r  e  V 
:(i)_ 


O) 


(15) 


The  set  of  equations  in  (15)  is  the  explicit  equivalent  of  the  equation  (14),  which  is  that  used  by  Backus 
and  Mulcahy. 

Backus  and  Mulcahy  then  define  the  stress  difference  —  S,j,  appearing  in  (15),  to  be  the  ‘stress 
glut’  r,j.  Next,  as  described  in  (v)  above,  they  linearize  the  equations  of  motion  everywhere  in  V,  the 
interior  of  the  planet,  so  that  in  (14)  they  in  effect  replace  pa,  by  p  d,2s,.  However,  as  we  sec  from  (15), 
this  requires  the  same  approximation  to  be  carried  out  in  the  nonlinear  region  V '.  That  is,  their  lineariza¬ 
tion  of  equation  (14)  throughout  the  planet  means  that  they  linearize  the  first  of  the  equatioas  in  (15) 


d 

r  <  -)  d 

,  ds, 

dt  1 

rh 

P  ”dT 

v.  J 

cannot  ordinarily 


be  approximated  by  p’3tV  in  such  a  region,  since  transport  terms,  which  may  be  large  compared  to  p’3t2s,’, 
arc  neglected  +.  However,  even  if  such  an  approximation  is  made  for  this  inert.  1  term  the  result  is: 


p  a2*,’  -  djS*1’  =  p  gj  +  a. 


S,<» 


reV 


(16) 


But  Backus  and  Mulcahy  use  the  camc  displacement  function  s,  everywhere  throughout  the  planet.  That 
is,  they  effectively  replace  the  inertial  term  p3,2s,  in  (16)  by  p^s/1’.  which  corresponds  to  the  analytic 
continuation  of  die  linear  (clcnic)  acceleration  field  into  the  nonlinear  failure  zone.  Therefore,  Backus 

^Transport  temis  are  those  arising  frum  the  last  term  tn  the  identity:  tl/clt  (pv)  =  c),  (pv)  +  vV(pv).  thus  when  the 
total  dertvati  .e  in  the  identity  is  approximated  hy  the  partial  derivative,  then  this  requires  neglect  of  icons  like  V'V(pv). 
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and  Mulcahy  use  a  fully  linearized  equation  in  the  failure  zone,  in  particular: 

P.a.V’  -  3jSf>  =  p.g,0’  *  3j(S,-  -  Sf)  ;  r  e  V '  (17) 

Thus,  these  authors  not  only  implicitly  linearize  the  inertial  term  in  equation  (15),  but  they  also 
replace  the  nonlinear  displacement  s,  with  the  analytic  continuation  of  the  linear  displacement  s,(1>  from 
outside  the  failure  zone  into  this  zone.  Further,  they  similarly  replace  p’  and  g'  by  p,  and  g/0,  again  by 
implicit  continuation.  Viewed  in  terms  of  the  validity  of  such  a  procedure,  even  as  the  crudest  of  approx¬ 
imations,  it  seems  clear  that  it  has  no  physical  or  logical  basis.  Certainly  if  fully  nonlinear  problems 
could  be  solved  by  introducing  a  source  term  which  incorporated  a  part  of  the  nonlinear  equation  (in  this 
case  the  nonlinear  stress  term),  with  the  remainder  of  the  equation  linearized  so  that  a  linear  Green’s  func¬ 
tion  solution  could  be  obtained,  then  we  would  have  at  our  disposal  a  method  that  could  solve  any  non¬ 
linear  problem.  There  is  no  question  that  this  is  not  the  case,  neither  in  general  nor  in  the  present  prob¬ 
lem. 

In  fact,  if  the  procedure  employed  by  Backus  and  Mulcahy  is  carried  to  its  correct  logical  conclu¬ 
sion,  then  the  differences  between  the  nonlinear  inertial  and  gravitational  terms  in  v'  and  the  analytic 
continuation  of  the  comparable  linear  terms  from  outside  the  failure  zone  into  V should  be  included  in 
the  source  factor.  That  is,  by  the  same  logic  used  to  introduce  the  stress  difference  term  (the  "stress  glut") 
in  the  equations  of  motion  in  V  ’,  one  can  form  from  (15)  (by  adding  and  subtracting  the  analytically  con¬ 
tinued  inertial  and  gravity  terms  p,^  s,(1)  and  Pig,(1))  the  equation: 


PA2  Si(l 


d  f-.l 

dT  H  ■ 


pAV 


pig,(1)  + 


gi(1>  +dj  S' 


This  is  perfectly  rigorous  and  no  approximations  arc  involved,  as  opposed  to  the  Backus-Mulcahy  result. 
Now  the  extra  factors  appearing,  in  addition  to  the  "stress  glut"  torn,  are  the  two  difference  terms  in 
brackets.  These  can,  following  the  reasoning  of  Backus  and  Mulcahy,  be  considered  as  (unknown) 
"source”  terms  and,  to  emphasize  this  interpretation,  can  be  written  in  combination  with  the  stress  glut 
term  as  a  "source"  of  the  linear  clastic  field  s,(1).  That  is,  the  previous  equation  can  be  rearranged  as: 
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Pi  d.2s,(1)  -  =  p,  g,( 


(i) . 


Pi  d&0)  -  —  (pi  v')  +  (p'g,' 


P.g.^  +  ^S'-S^)  ;  teV 


Further,  one  can  combine  this  result  with  the  second  equation  in  (15)  to  obtain  a  result  applicable  in  "all 
space",  that  is  for  r  e  V  ©V).  In  particular,  one  has  the  formal  result: 


Pi  3t2  s,(i)  -  d]  S<’>  =  p,g,(1)  +  y,  :  r  e  V  ©  V, 


where 


Y,  =  - 


Pi^s'0-  ^(pvj 


p  g,  -  pig; 


o) 


+  9J(S0-S<1V  ;reV 


Y.  =  0  ;  r  e  V, 

This  is  essentially  the  same  kind  of  result  as  that  obtained  by  Backus  and  Mulcahy,  but  includes  the  two 
extra  factors  in  the  source  term  y,  that  were  neglected  by  them. 

At  this  point  one  might  be  tempted  to  conclude,  since  this  is  rigorous  with  nothing  neglected,  that 
the  result  constitutes  a  very  similar  but  somewhat  more  "precise"  representation  of  the  "equivalent 
source"  term.  However,  while  it  is  true  that  the  equation  is,  in  fact,  the  correct  consequence  of  the  pro¬ 
cedure  used  by  Backus  and  Mulcahy,  it  by  no  means  provides  a  representation  of  the  equivalent  source 
due  to  spontaneous  failure  (an  earthquake)  or  any  other  kind  of  source.  This  is  easy  to  see,  since  it  is  only 
necessary  to  note  from  the  first  of  the  equations  in  (15)  that: 

|-(p'v')-ajS,;  =  p  g‘  ;  r  e  V  ’ 

which  only  expresses  conservation  of  momentum  in  V  .  But  if  this  equality  must  be  satisfied,  which  of 
course  it  must,  then  several  of  the  terms  in  the  expression  for  y,  must  cancel.  In  particular,  using  this 
equation  we  have  that  y  must  always  (identically)  reduce  to: 

Y,  =  jp|3,2s,(1)  -  djS.j0  -  p,g,(,>|  ;  r  e  V  ' 

Y,  ~ U  ;reV, 

If  this  is  now  re-inserted  into  the  "equation  of  motion  in  all  space",  we  have: 
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Pi^V0  -  djSf  =  Pigi0)  +  jpiS/s/0  -  3jS<1)  -  p,g,(1)j  ;  r  e  V ' 
p,ai2s1(1)-3jS/1)  =  Pig/°  ;rev, 

While  these  equations  are  exact,  the  first  is  nothing  but  a  trivial  identity  and  the  second  is  the  regu¬ 
lar  equation  of  motion  in  the  elastic  medium  outside  the  source  region.  Thus,  the  consequence  of  the 
approach  used  by  Backus  and  Mulcahy,  when  carried  out  properly,  is  a  triviality.  However,  since  they 
actually  neglect  terms  of  the  same  order  as  the  "stress  glut"  term  they  retain,  they  do  not  recognize  the  cir¬ 
cular  nature  of  their  procedure.  In  any  case,  it  is  clearly  safe  to  conclude  that  the  formal  mathematical 
basis  for  their  representation  is  erroneous. 


IV.  Consequences  of  the  Stress  Glut  Phenomenological  Representation 

Aside  from  the  previous  observations  concerning  the  foundations  of  the  stress  glut  representation, 
there  arc  several  consequences  of  this  formulation  that  indicate  its  inappropriate  nature. 

As  a  first  example,  consider  the  limiting  form  of  the  stress  glut  function  as  time  tends  to  infinity 
(the  static  limit).  In  this  case,  with  y,v  =  3)riJ: 


lim  y*  =  lim 

t~>oo  l— 


-  lim  [W]  =| 

L  J 

pg>-plg. 
0;  r  e  V, 


(F)  . 


;  r  e  V 


where  p\  etc.,  represent  final  static  values  of  density,  etc.  Here  we  have  simply  used  the  fact  that  the 
static  stresses  must  satisfy  equilibrium  equations.  Therefore  we  observe  that,  like  Gilbert's  representa¬ 
tion,  this  source  term  only  represents  gravity  field  changes  in  the  long  time,  or  low  frequency,  limit. 
However,  while  Gilbert’s  result  represented  earthquake  induced  gravity  changes  in  the  entire  planet,  the 
result  given  by  Backus  and  Mulcahy  represents,  in  the  low  frequency  limit,  only  the  changes  within  the 
failure  zone.  But,  as  previously  observed,  the  total  of  all  gravity  change  effects  arc  at  least  several  orders 
smaller  than  tectonic  effects  and  the  effect  of  those  within  the  failure  zone  alone  would  be  considerably 
less  than  those  to  be  obtained  from  Gilbert’s  result. 
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More  explicitly,  if  we  evaluate  the  displacement  field  due  to  the  "stress  glut",  as  given  by  Backus 
and  Mulcahy  in  equation  (1 1),  taking  the  stress  functions  involved  in  the  stress  glut  tensor  to  be  at  least 
approximately  separable  in  space  and  time,  then  we  get: 


sk(1)(x,o=£ 

v=l 


J(pg,j-p,g/1))u;dv 


J*  sin  u)v(t  -  T)  f  (t)  dx 
o 


where  f(t)  denotes  a  time  function  characterizing  the  temporal  variation  cf  both  the  stresses  s’  and  S,j  in 
the  definition  of  the  stress  glut.  Thus,  no  matter  what  the  form  of  the  lime  variation  of  the  stress  glut  fac¬ 
tor,  the  magnitude  of  the  predicted  displacement  field  is  determined  by  the  inner  product  of  the  gravity 
field  changes  with  the  elastic  eigenfunctions  for  the  earth  taken  over  the  failure  zone.  As  with  Gilbert’s 
result,  with  wnich  it  differs  only  in  that  the  volume  integration  is  over  the  failure  zone  rather  than  the 
entire  volume  of  the  earth,  the  predicted  displacement  field  from  this  representation  will  be  many  orders 
of  magnitude  less  than  that  actually  associated  with  earthquakes. 

The  use  of  the  stress  glut  formulation  in  the  context  of  inclusion  theory  (c.g.  Eshelby,  1957),  affords 
another  example  of  a  clearly  incorrect  prediction  for  the  clastic  wave  radiation.  Here  the  problem,  as  con¬ 
sidered  by  Backus  and  Mulcahy,  involves  the  creation  of  a  volume  inclusion  within  a  stressed  medium 
where  the  material  within  the  inclusion  is  viewed  as  transforming  to  a  new  physical  state,  such  as  from  a 
solid  to  a  fluid  +.  A  transformation  of  this  sort  must  result  in  the  radiation  of  elastic  waves,  since  the 
stress  outside  the  inclusion  must  dynamically  adjust  to  the  presence  of  the  new  material.  In  treating  this 
problem  Backus  and  Mulcahy  use  the  expression 

fy  =  E|ji,lckl 

for  the  stress  glut,  where  Eljkl  is  taken  to  be  the  clastic  tensor  throughout  the  inclusion  after  transformation 
and  c£  is  the  so-called  "stress-free  strain",  which  is  precisely  the  strain  that  would  occur  in  the  inclusion  if 
it  were  allowed  to  undergo  the  transformation  unconstrained  by  the  surrounding  matrix  {i.c.,  when  the 

+  The  lime  variation  for  this  transformation  is  unspecified  in  the  problem 
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inclusion  boundary  is  a  free  surface). 

However,  for  the  transformations  described  by  the  Eshelby  inclusion  theory,  must  be  a  constant 
(Eshclby,  1957).  Further  c£  is  by  definition  an  equilibrium  field.  Consequently,  the  equivalent  volume 
force  vector  associated  with  this  stress  glut  vanishes.  That  is: 


since  the  clastic  tensor  is  constant  in  the  inclusion  and  e£  is  an  equilibrium  strain  field  by  construction. 
Therefore  this  stress  glut  representation  gives  a  null  result  for  this  problem  and  this  obviously  cannot  be 
correct. 

It  is  also  easy  to  see  that  a  stress  glut  formulation  does  not  describe  the  physical  manifestations  of 
phase  transitions.  In  particular,  consider  the  transition  to  involve  sudden  melting  to  form  an  inclusion 
within  a  solid  matrix  that  is  initially  under  a  pure  shear  prestress.  Suppose  the  melting  is  such  that  an 
ideal  fluid  is  produced,  so  that  the  shear  modulus  is  zero.  Then  Eijkl  reduces  to  a  form  such  that 

r,j  ~ 

where  X  is  the  Lam6  constant  for  the  fluid  and  is  the  volume  dilation  of  the  fluid  when  the  boundary' 
stresses  due  to  the  surrounding  matrix  arc  removed.  From  the  form  of  the  stress  glut  it  is  evident  that  it 
corresponds  to  a  pure  pressure,  with  no  shear  stress  components.  Thus,  we  can  conclude  that  this  source 
can  never  produce  a  correct  result,  since  not  only  docs  the  divergence  of  the  quantity  r,  vanish  but  the 
only  radiation  field  that  could  possibly  be  produced  by  such  a  source  term  would  be  a  pure  comprcssional 
wave  field  since  only  is  present  and  wc  know  that  shear  waves  must  be  produced  because  the  initial 
field  in  the  surrounding  matrix  was  pure  shear  and  it  must  change  drastically  after  the  transformation. 
(Such  a  change  can  only  be  accomplished  by  the  radiation  of  shear  waves,  which  result  in  a  relaxation  of 
the  shear  stresses  in  the  vicinity  of  the  fluid  inclusion.)  Therefore,  aside  from  the  earlier  fundamental  con¬ 
siderations,  even  a  rather  cursory  examination  of  these  special  cases  shows  that  the  stress  glut  source 
representation  fails  to  describe  the  ba.->ic  properties  of  the  clastic  wave  radiation  to  be  expected. 
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V.  Dynamical  Representation  of  Spontaneous  Failure  in  a  Stressed  Medium:  The  "Time  Depen¬ 
dent  Equilibrium"  Method 

In  order  to  display  the  proper  "equivalent  force"  term  for  a  spontaneous  seismic  source  like  an  earth¬ 
quake,  which  derives  energy  for  its  formation  from  the  medium  itself  and  simultaneously  releases  stored 
energy  from  the  medium,  as  evidenced  by  the  seismic  waves  radiated  from  the  event,  it  is  necessary  to 
carefully  consider  the  nature  of  equilibrium  in  the  medium  during  formation  of  a  failure  zone.  In  doing 
so  we  observe  from  the  onset  that  the  conditions  for  equilibrium  arc  continually  changing  with  time  as  the 
failure  z.onc  forms  and  expands.  This  follows  from  the  force  balance  requirements  themselves  and,  in 
particular,  from  the  fact  that  equilibrium  depends  upon  continuity  conditions  at  all  the  boundaries  of  the 
medium.  Therefore  since  failure  produces  a  boundary  within  the  medium  that  changes  dimensions  with 
time,  then  the  equilibrium  stress  must  also  change  parametrically  with  time  because  of  the  changing 
boundary  conditions  applying  to  the  internal  stress  fields. 

More  quantitatively  these  statements  regarding  equilibrium  changes  are  expressed  by  the  rela¬ 
tions  +: 

' 

di4e)  =  -  P(e)gke) ;  r,  e  V,(0 

(18) 

Q/4eVll  =  0;r,e  av^g 

where  we  use  a  superscript  (e)  to  explicitly  label  the  variables  as  equilibrium  fields.  Further,  we  choose 
to  explicitly  label  the  coordinates  as  r,  and  t,,  denoting  "source  spatial  coordinates”  and  a  "source  lime 
coordinate",  respectively.  These  equations  arc  entirely  analogous  to  those  of  equation  (2),  except  that  we 
have  taken  explicit  note  of  the  fact  that  the  volume  external  to  the  failure  zone,  V,  as  shown  in  Figure  (I), 
and  the  surface  of  this  volume,  dv,,  are  both  parametrically  dependent  on  the  time  variable  t*.  Thus  the 
equations  in  (18)  express  the  conditions  for  equilibrium  in  the  clastic  region  and  since  V , and  dV,(i,) 
depend  on  time,  then  the  solution  of  these  equations  will  yield  stress,  density  and  gravity  fields  that  arc 

+As  before,  ihe  double  bracket  notation  appearing  in  (18)  is  used  to  denote  the  jump  in  the  enclosed  quantity  across  the 
boundary. 
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functions  of  both  r,  and  t,.  This  dependence  will  occasionally  be  explicitly  displayed,  for  example  by 
writing  the  stress  x^  as  xic)(r..  t,)  and  similarly  for  the  other  field  variables.  Likewise,  the  equations 
specifying  equilibrium  inside  the  failure  zone  are 

a/1)x«  =  -p(V;rie  V'(o] 

I(1)tie)nj  =  0;r,€  3V'(g  |  ^ 

where  V'  is  the  failure  zone  volume  and  3V'  its  surface,  as  shown  Figure  (1).  This  equation  expresses  the 
same  equilibrium  relationships  as  were  expressed  by  equation  (3),  but  here  again  we  have  been  more 
explicit  in  indicating  the  parametric  dependence  of  the  fields  on  time  through  the  linkage  with  the  time 
dependent  volume  and  surface  boundary  of  the  failure  zone. 

Of  course  it  is  not  usual  to  speak  of  an  equilibrium  field  as  depending  on  time  and  it  appears  at  first 
sight  to  be  a  contradiction  in  terms.  However  there  is  really  no  contradiction,  since  the  dependence  is 
parametric;  that  is  these  fields  are  the  equilibrium  fields  that  would  exist  if  the  failure  zone  boundary  were 
of  the  size  and  shape  specified  at  a  particular  time  t,.  Thus,  we  are  merely  asking  what  the  equilibrium 
fields  would  be  if  the  failure  zone  were  "frozen",  for  all  time,  at  a  size  and  shape  appropriate  to  a  particu¬ 
lar  time  t,.  Since  the  size  and  shape  of  the  failure  zone  is  a  function  of  the  time  variable  t,  and  since  the 
equilibrium  fields  (stress  or  displacement)  arc  a  function  of  coordinates  of  the  surface  at  time  t,,  then  the 
equilibrium  fields  will  be  implicit  functions  of  t,.  Further,  we  know  that  the  spatial  dependence  of  an 
equilibrium  displacement  field  is  always  a  bi-harmonic  function  with  coefficients  oi  the  bi-harmonic 
scries  depending  on  the  surface  coordinates  of  the  failure  zone  inclusion.  Thus  we  also  know,  a-priori, 
that  the  parametric  dependence  on  t,  will  appear  in  the  coefficients  of  a  bi-harmonic  series  expansion  for 
the  equilibrium  displacement. 

In  the  dynamic  description  of  the  phenomena  associated  with  spontaneous  failure  in  a  stressed 
medium  the  time  dependence  of  the  equilibrium  field  is  a  hidden  (or  internal)  variable.  That  is,  the  equa¬ 
tions  of  motion  in  the  linear  region  outside  the  failure  zone  are 
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ajp  a,uk)  -  akxlk = p  gk ;  k  =  i ,  2, 3 


(20) 


when  written  with  the  source  coordinates  as  independent  variables.  However,  the  various  dependent  vari¬ 
ables,  such  as  the  displacement  field  uk(r„  t,),  arc  total  fields;  that  is  they  include  both  the  equilibrium 
field  and  a  purely  dynamic  field  as  components.  In  particular,  they  arc  of  the  form 

uk  =  uk(e)  +  uk(d)  (21) 

where  the  "dynamic  field",  ukd),  can  be  defined  simply  as  the  difference  between  the  field  uk  appearing  in 
(20)  and  the  field  u^e)  defined  by  (18).  Clearly,  (20)  docs  not  display  the  equilibrium  field  uk(e)  explicitly 
and  so  it,  along  with  the  other  equilibrium  fields  connected  with  this  equation,  can  be  termed  hidden 
internal  field  variables. 

Ordinarily  the  equilibrium  fields  are  of  no  consequence  in  the  solution  of  (20),  since  in  most  elasto- 
dynamic  problems  the  equilibrium  displacement  is  not  a  function  of  time.  Thus,  in  the  case  when  the 
equilibrium  is  constant,  substitution  of  the  complete  form  for  uk  into  (20)  produces  djp  duuk(d))  for  the 
inertial  term  and  substitution  of  the  comparable  expressions  for  %  and  gk,  involving  the  sum  of  their 
equilibrium  and  dynamic  components,  results  in  cancellation  of  the  resulting  equilibrium  terms  by  virtue 
of  the  equations  of  equilibrium,  as  given  by  equations  of  the  form  of  (18).  Therefore  in  the  case  of  time 
independent  equilibrium,  equation  (20)  is  an  equation  for  the  dynamic  fields  and  states  that: 

a,(P  auuk(d))  -  a,4d>  =  P  gk(d) ;  k  =  1 , 2, 3 

However,  when  the  equilibrium  fields  arc  lime  dependent,  as  is  the  situation  for  a  spontaneous 
failure  process,  then  equation  (20)  gives  a  different  result.  In  this  case  wc  have 

uk  =  uke)(r„  O  +  uk(d)(rs,  Q 

hk=4')(r..O  +  4d)(r,,0-  (22) 

gk  =  gkC)(rs.  O  +  gkd)(rs>  O 

where  the  equilibrium  fields  all  depend  on  t,  by  virtue  of  equations  (18)  and  (19).  Now,  use  of  these 
forms  in  (20)  gives: 
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du(p  awuk(d))  +  a^p  auukw)  -  a,^  -  a^  =  p  gk(e) + p  g^d) 


Here  the  inertial  term  involving  u^e)  does  not  necessarily  vanish,  since  there  will  be  variations  in  uk<e) 
because  of  the  time  dependent  boundary  changes.  However,  as  previously  noted  in  equation  (18): 

di4e)=-pgke) 

and  this  relation  holds  in  all  situations;  that  is  whether  and  are  parametrically  time  dependent  or 
not.  Therefore  we  have, 

3U(P  3uukd))  -  ^ikd)  =  -  ^U(P  K^c)) +  P  gkd) ;  r,  e  V,(0 

The  equilibrium  field  in  the  linear  zone  outside  the  failure  volume  can,  ho  ■  ever,  be  expressed  as  the 
sum  of  the  changes  in  the  equilibrium  state  due  to  the  introduction  of  the  failure  zone,  which  can  be 
denoted  by  u(e)'  (r,  ,  4),  plus  the  (time  independent)  field  u®  (rj  describing  the  initial  equilibrium  state. 
Thus  we  can  always  express  u(e)  as 

u(e) (r,,t,)  =  u<e>'(rt  ,  Q  +  u®(r,)  ;  r.eV, 

as  a  consequence  of  linearity.  Since  time  derivatives  of  u®  vanish,  then  the  equation  of  motion  becomes 

djp  auuk(d))  -  a,  x£d)  =  -  a,,  (p  au  uk(e)')  +  pgk(d)  ;  r,  e  V,(t,)  (23) 

Thus,  we  have  that  the  term  -  at(p  auuk(e)  )  is  the  source  term  (or  "equivalent  body  force  term")  for  the 
dynamic  field  generated  by  a  localized  failure  process  within  a  prestressed  solid.  Clearly  this  source  term 
is  distributed  throughout  the  region  V„  since  it  involves  time  rate  changes  of  the  equilibrium  displace¬ 
ment  defined  throughout  this  region.  Thus,  the  equivalent  body  force  term  for  an  earthquake  is  not  the 
divergence  of  a  stress  drop,  but  is  an  "inertia!  force"  produced  by  changes  in  the  elastic  equilibrium  dis¬ 
placement  field  that  occur  because  of  the  creation  and  growth  of  the  failure  zone  in  the  prestressed  solid. 
More  specifically,  if  the  failure  zone  growth  occupies  a  time  interval  (0,  Tq),  then  neglecting  the  small 
effects  of  density  changes  associated  with  the  relaxation  of  stresses: 

O  ;  t,  <  0 

-  a,(p  auuk<c)')  =  -  p  aw2uk(e)' = •  F(r„  g  *  o ;  o  <  t,  <  t0 

0 ;  t,  >  t0 
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where  the  "force"  function  F(rs,  ts)  is  a  direct  function  of  the  failure  volume  lime  history  of  growth.  This 
function  must  also  depend  on  the  rheology  of  the  material  within  the  failure  zone  as  well  as  the  elastic 
properties  of  the  medium  surrounding  it,  by  virtue  of  the  boundary  conditions  at  the  failure  surface. 

The  analytical  form  of  this  source  term  is  prescribed  by  the  boundary  value  problem  of  (18).  There¬ 
fore,  since  the  fields  prescribed  by  (18)  are  equilibrium  fields,  it  follows  that  the  displacement  field  u(c)  is 
a  biharmonic  function  of  the  spatial  coordinates  (c.g.  Landau  and  Lifshitz,  1951),  as  noted  earlier. 
Further,  the  field  u(I)  and  u(e)  arc  similarly  biharmonic.  Thus,  since  the  coefficients  of  the  biharmonic 
scries  arc  fixed  by  the  shape  of  the  failure  zone  inclusion  as  well  as  by  the  value  of  the  initial  stress  field 
and  the  rheology  of  the  material  (both  within  and  outside  the  failure  zone),  then  the  parametric  time 
dependence  of  u(e)  will  be  contained  within  these  series  coefficients.  Consequently,  even  without 
knowledge  of  the  prestress  field,  the  rate  of  failure  expansion,  the  material  properties  and  the  shape  of  the 
failure  zone,  we  know  a  good  deal  about  the  analytical  form  of  this  equivalent  force  term. 

Before  making  use  of  this  formulation  for  the  dynamics  of  failure  phenomena  (and  our  knowledge 
of  the  biharmonic  form  of  u(c)  )  it  is  useful  to  point  out  that  this  formulation  is  a  particular  example  of  a 
more  general  class  of  problems  involving  time  dependent  internal  equilibrium  changes.  In  this  regard,  it 
has  been  shown  by  Dilts  (1985)  that  the  decomposition  of  stress  and  displacement  fields  into  dynamic  and 
equilibrium  components,  so  as  to  display  the  internal  equilibrium  field  variable,  allows  the  permanent 
deformations  due  to  microcracking  or  lattice  dislocation  creation  or  movement  (which  results  in  macros¬ 
copic  plastic  behavior)  to  be  treated  directly  in  terms  or  equilibrium  field  changes  in  the  equations  of 
motion  for  the  material.  As  with  the  present  problem,  time  dependent  equilibrium  changes  appear  as 
dynamical  terms  in  the  equations  of  motion.  However,  as  shown  by  Dills,  since  the  particular  equilibrium 
field  changes  associated  with  microscopic  disordering  within  the  solid  changes  arc  generally  dependent 
on  the  dynamical  displacement  u(d>  and  its  derivatives,  then  the  result  is  the  occurrence  of  non-linear 
terms  in  the  equations  of  motion. 
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On  the  other  hand,  we  observe  that  the  result  expressed  by  (23)  is  a  linear  equation.  This  is  because 
we  can  partition  the  spatial  region  around  a  single  failure  zone  into  zones  of  linear  and  nonlinear  material 
behavior,  where  the  changes  in  equilibrium  within  the  linear  (clastic)  zone  only  occur  because  of  the 
requirement  that  boundary  conditions  must  be  satisfied  on  the  time  varying  failure  boundary.  Therefore, 
the  equilibrium  field,  u(e) ,  at  a  point  in  V[  docs  not  depend  on  the  dynamic  displacement  or  its  derivatives 
at  that  point,  but  only  on  geometric  and  material  characteristics  at  the  failure  boundary.  As  a  conse¬ 
quence  the  equilibrium  term  in  (23)  is  linear.  However,  the  behavior  of  the  material  within  the  failure 
zone  would  be  described  by  internal  equilibrium  field  changes  that  depend  on  the  local  displacement  field 
and  its  dcratives  due  to  the  intense  microscopic  disordering  produced  by  any  failure  process.  Conse¬ 
quently  the  dynamics  of  the  material  motion  in  the  failure  zone  would  be  described  by  an  equation  of 
motion  with  dynamic  terms  again  involving  changes  in  internal  equilibrium,  but  now  with  these  terms 
depending  on  the  local  dynamic  displacement  and  its  derivatives,  so  that  these  terms  would  be  nonlinear 
in  the  dynamic  field  variable. 

Wc  note  further  that  both  linear  and  non-linear  changes  in  equilibrium  can  result  in  either  absorp¬ 
tion  of  energy  from  the  dynamic  field  or  release  of  energy  to  it.  That  is,  ’’dynamic  equilibrium"  terms  in 
the  equations  of  motion  can  correspond  to  either  "sources"  or  "sinks"  of  energy  for  the  dynamic  field. 
Further,  this  behavior  may  not  only  be  spatially  dependent  but  also  time  dependent,  so  that  energy  may  be 
absorbed  or  stored  at  a  point  and  later  partially  or  totally  released  at  that  same  point.  Clearly,  because  of 
the  irreversible  nature  of  the  nonlinear  changes  within  a  failure  zone,  we  expect  absorption  of  energy  in 
this  region  while  in  the  linear  zone  wc  expect  a  net  reduction  of  stored  energy  in  the  equilibrium  field, 
with  the  stored  energy  changes  being  a  consequence  of  relaxation  of  the  equilibrium  tectonic  prestress 
field  around  the  failure  zone.  Clearly  then,  the  contribution  of  the  equilibrium  field  changes  as  "sources" 
for  the  dynamic  field  can  be  very  complex  and  can  involve  temporary  increases  in  stored  energy,  followed 
by  decreases,  since  the  strain  changes  arc  reversible  in  the  linear  zone  and  the  growing  failure  zone  can 
require  increased  strain  levels  locally  followed  by  decreases  as  the  position  of  the  point  relative  to  the 
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failure  boundary  changes.  Such  behavior  will  be  dictated,  on  a  spatial  and  temporal  basis,  by  die  analyti¬ 
cal  fonm  of  u(e)  in  equation  (23),  with  the  time  dependent  biharmonic  coefficients  in  the  expansion  for  u(c) 
reflecting  the  changing  position  of  any  point  in  the  medium  relative  to  the  failure  boundary  as  it  grows. 


In  addition  to  the  rather  complex  linear  behavior  to  be  expected  from  the  volume  source  term  in 
equation  (23),  it  is  also  to  be  expected  that  dynamical  complications  will  arise  at  the  boundary  of  the 
failure  zone.  Here  we  expect,  on  physical  grounds,  that  energy  will  be  transferred  from  the  linear  zone  to 
the  nonlinear  zone  where  it  will  be  absorbed  by  one  irreversible  process  or  another.  In  particular,  it  can 
be  expected  that  the  energy  required  to  cause  failure  will  be,  at  least  in  part,  extracted  from  the  dynamical 
field  at  the  boundary.  Such  phenomena  must  therefore  be  reflected  in  the  boundary  conditions  connecting 
the  linear  and  non-linear  zones  at  the  failure  boundary.  Thus  such  boundary  conditions,  along  with  (23), 
serve  to  define  the  complete  dynamical  boundary  value  problem  for  the  displacement  field  in  the  linear 
region  V,. 

The  boundary  conditions  that  arc  appropriate  are  those  that  express  conservation  of  mass,  momen¬ 
tum  and  energy  across  such  a  "singular"  surface,  where  the  possibility  of  irreversible  energy  absorption 
along  a  moving  boundary  is  explicitly  addressed.  This  physical  situation  has  been  considered  by 
Archambcau  and  Minster  (1978)  and  they  have  shown  that  the  appropriate  boundary  conditions  on  the 
failure  surface,  3V',  are: 

Ipv.*nJ  =0 

E(pvkv*  -  Tkl)n,]j  =  0  •  (24) 

([(p£  v  *  -  vkxkj  +  q,)n,]J  =  0 

- 

where  the  double  bracket  denotes  the  change  in  the  quantity  enclosed  across  the  failure  surface.  Thus,  for 
example: 

ffvj  *  v.(dV,V  v,(dVo0 

where  3V,'  and  dv0'  denote  the  surface  av '  approached  from  inside  and  outside  the  failure  zone. 
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In  these  relations  v  denotes  the  particle  velocity  —  ,  n  is  the  normal  to  the  failure  zone  boundary 


V.  J 

and  q  the  heat  flux  vector.  Also 

v*s  v  -  U 

where  U  is  the  velocity  of  the  failure  surface  boundary  and  E  is  the  total  energy  density  for  the  material, 
defined  by 

£  s  T  +  vkvk/2  + 

with  r  the  internal  energy  and  <j>  the  gravitation  energy.  The  quadratic  vkvk/2  corresponds  to  the  kinetic 
energy  contribution. 

All  the  field  variables  appearing  in  (24)  are,  of  course,  functions  of  the  coordinates  and  time,  includ¬ 
ing  the  surface  normal  n  and  U,  the  rate  of  change  of  the  failure  boundary  with  time.  The  first  of  the 
equations  in  (24)  expresses  conservation  of  mass,  the  second  conservation  of  momentum  and  the  last  con¬ 
servation  of  energy.  All  the  field  variables  include  equilibrium  as  well  as  dynamic  components.  Thus, 
for  example,  the  particle  velocities  have  the  explicit  form: 

v  =  v<d>  +  v(e) 

As  was  shown  by  Archambcau  and  Minster,  the  boundary  conditions  can  be  reduced  to  simpler 
form  lor  the  physical  situation  prevailing  during  rapid  failure.  In  particular,  when  the  density  changes 
(and  the  gravitational  change  <J>)  upon  failure  are  small  and  the  rupture  rate  is  high  (that  is  equal  to  a 
significant  fraction  of  the  local  shear  velocity  and  much  larger  than  the  particle  velocity  v)  then,  when  the 
heat  flux  changes  across  the  failure  surface  are  neglected  relative  to  the  larger  mechanical  flux  term 

l[vkTk,]],  it  follows  that  the  equations  in  (24)  reduce  to: 

pUR[[vk]  =  -  [[tJ  (25) 

where  tk  -  tk,n,  arc  the  components  of  the  surface  tractions  on  the  failure  surface.  Equation  (25)  combines 
the  conservation  of  mass  and  momentum  equations  from  (24).  Further,  under  rapid  failure  conditions  the 


energy  equation  in  (24)  reduces  to: 


IvJ  =  2p2Uk3L 


(26) 


where  the  change  in  internal  energy,  m  ,  is  identified  as  a  material  property  (ie.  a  "latent  heal"  of  transi¬ 
tion)  and  is  denoted  as  L.  Here  the  "rupture  rate",  UR,  is  defined  as 

UR  3  (U  -  v)  n 

and  expresses  the  excess  of  the  expansion  rate  of  the  failure  surface  over  the  particle  velocity  in  the  direc¬ 
tion  normal  to  the  surface. 

Equation  (26)  expresses  the  fact  that  an  energy  barrier  to  failure  exists,  and  this  barrier  is 

represented  by  the  internal  energy  change  [[r]]  (or  L)  required  for  failure  to  occur  in  the  material. 
Further,  the  equation  shows  that  the  rupture  rate  is  proportional  to  the  magnitude  of  the  traction  change 
(the  "stress  drop")  divided  by  square  root  of  this  energy  factor,  that  is: 


Smee  p  and  L  arc  material  constants,  then  this  says  that  the  magnitude  of  the  traction  or  stress  "drop", 
1/2 

I  At  |  s  HvkU  .  is  directly  proportional  to  the  rupture  rate  and  vice-versa.  This  is  a  dynamical  condition 
that  will  reflect  itself  in  the  proper  solutions  for  radiation  fields  from  failure  processes  in  different 
material  types. 

The  two  dynamical  boundary  conditions  on  the  failure  surface  can  be  combined  so  as  to  eliminate 
UR  and  produce  a  single  boundary  condition  involving  only  the  dynamic  field  variables  and  the  intrinsic 

internal  energy  change,  represented  by  the  (material)  parameter  L.  That  is  from  (25)  and  (26)  we  have: 

1/2 

M«-[2L/ffvjj  M  (27) 

assuming  that  none  of  the  traction  jumps  vanish  (in  which  ease  the  conservation  relations  are  trivially 
satisfied.) 
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Clearly  this  relation  shows  that  the  jump,  or  discontinuity,  in  the  particle  velocity  vector  across  the 
failure  boundary  is  in  the  opposite  direction  as  the  jump  in  the  traction  vector  (and  vice-versa)  and  is,  as 
well,  proportional  to  the  square  root  of  the  material  parameter  L.  However,  the  relationship  between  par¬ 
ticle  velocity  and  traction  jumps  is  non-linear,  due  to  the  dependence  on  the  traction  ratio  on  the  right  side 
of  (27).  Thus,  while  (25)  expresses  a  similar  working  relationship  between  the  traction  and  velocity 
jumps  across  the  failure  surface,  the  true  non-linear  relationship  between  these  field  variables  is  hidden  in 
the  factor  involving  the  rupture  rate.  In  particular,  since  the  rupture  rate  is  required  to  be  proportional  to 
the  stress  drop  factor  |  At  |  in  order  that  energy  be  conserved,  it  follows  that  (25)  represents  a  non-linear 
relationship  between  the  traction  and  particle  velocity  field  variables. 

It  is  worth  noting  that  the  dynamics  of  the  failure  process  requires  that  both  the  particle  velocity  and 
the  tractions  be  discontinuous  across  the  failure  boundary.  This  is  not  the  usual  assumption  made  in  the 
construction  of  pure  kincmatical  models  of  failure,  such  as  in  a  dislocation  model  of  an  earthquake  where 
the  tractions  are  assumed  to  be  continuous  while  components  of  the  particle  velocity  (and  displacement) 
arc  assumed  discontinuous.  Neither  is  it  the  case  for  the  variety  of  so  called  stress  pulse  models,  where  a 
stress  drop  on  the  failure  surface  is  used  to  "drive"  the  surrounding  (elastic)  medium  and  thereby  produce 
a  simulated  earthquake  radiation  field.  Obviously  when  (current)  dislocation  or  stress  drop  models  are 
used  to  generate  a  moment  representation,  for  purposes  of  inverting  observational  data  from  earthquakes, 
the  result  will  be  non-physical  in  the  sense  that  the  boundary  conditions  assumed  are  not  compatible  with 
the  dynamical  constraints  imposed  by  energy  and  momentum  conservation.  Therefore  "displacement 
offsets"  or  "stress  drops"  inferred  from  observations  ob'aincd  through  the  use  of  such  formulations  do  not 
actually  represent  these  physical  quantities.  As  is  shown  in  a  later  section  (VIII),  the  required  form  of  the 
Green’s  function  surface  integral  representation  describing  the  seismic  radiation  involves  both  traction 
and  displacement  changes  on  the  failure  surface.  Therefore  a  combination  of  a  "stress  pulse"  nd  "dislo¬ 
cation"  sources  is  generally  required  for  compatibility  with  the  dynamical  constraints  involved  in  spon¬ 
taneous  failure. 
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VI.  The  Green’s  Function  Representation  of  Dynamic  Fields  Generated  by  Time  Dependent  Equili¬ 
brium  Changes 

A  compact  representation  of  the  dynamic  field  u(d)  appearing  in  the  equations  of  motion  (23),  which 
also  incorporates  the  effects  of  the  boundary  conditions  in  (25)  and  (26),  is  obtained  by  introducing  the 
usual  Green’s  function  for  the  clastic  medium 

Specifically,  consider  the  Green’s  function  Gkm  (r,  t;  r^t,)  satisfying  the  inhomogeneous  clastic  wave 
equation: 

di,(p3t,Gkm)-diG£  =  4jt  8£"  8(r  -  r,)8(t  -  q)  =  4nA™  (28) 

where  G™  is  a  two  point  tensor  obeying  the  usual  causality  relations  (see  Archambeau  and  Minster,  1978, 

for  details)  and  where  r,  and  t,  are  source  coordinates  and  time,  while  r  and  t  arc  receiver  coordinates  and 
time.  Here  8^  is  the  Kronccker  delta  (which  is  zero  if  k*m  and  unity  if  k  =  m),  while  8(r  -  r,)  and  8(t  - 1,) 
are  Dirac  delta  functions.  All  partial  derivatives  in  (28)  are  to  be  taken  with  respect  to  the  source  coordi¬ 
nates.  Further,  G£  denotes  the  stress  tensor  associated  with  Gkm.  That  is: 

G%=  c1Jki  a.c; 

where  Cijk)  is  the  fourth  order  elastic  tensor. 

Wc  can  take  Gkm  to  satisfy  boundary  conditions  along  the  exterior  surface  of  the  medium  and  along 
interior  surfaces  of  material  discontinuity,  excluding  the  (growing)  failure  surface ,  since  in  the  cases  of 
most  interest  explicit  representations  of  G™  in  layered  media  are  known.  That  is,  G™  may  be  taken  to 
satisfy  all  "normal"  boundary  conditions  within  flat  or  spherically  layered  media  except  those  "non¬ 
standard"  conditions,  given  by  (25)  and  (26),  which  apply  to  the  failure  boundary.  An  explicit  form  for 
Gkm  which  satisfies  the  normal  material  boundary  conditions  involving  continuity  of  tractions  and  the  nor¬ 
mal  component  of  the  particle  velocity  (or  all  particle  displacement  components  if  the  material  boundary 
is  a  "welded"  solid-solid  interface)  is,  in  the  frequency  domain: 


G,m(r,  co;r,,co,)  =  4n]£ 
k(u) 


\|/m(r.co)  n/|(r,.0),) 
(co2  -  co2)N(k,  co) 


(2b) 
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where  ym  are  the  component  displacement  eigenfunctions  for  the  medium  ( without  the  failure  zone)  and 
k((o)  denotes  the  infinite  set  of  eigenvalues  for  the  medium.  (Here  y,  is  the  complex  conjugate  of  y,  and  N 
(k,  o) )  is  a  normalization  constant.)  Since  the  eigenfunctions  y  satisfy  the  usual  material  boundary  condi¬ 
tions  by  definition,  then  G/”  also  does.  For  the  eigenfunctions  in  layered  spheres  and  half  spaces  see,  for 
example,  Ben  Menahem  and  Singh  (1972,  1981)  or  Harvey  (1981).  (For  a  layered  half  space,  however, 
the  sum  over  the  wave  number  k(co)  in  the  eigenfunction  expansion  for  Gjm  may  be  replaced  in  pan,  or 
totally,  by  an  integral,  so  the  summation  in  (29)  should  be  interpreted  as  a  generalized  sum.)  These  simple 
"layered  media"  cases  are  those  of  most  interest  in  geophysics  and  in  any  case  are  sufficient  for  a  study  of 
source  radiation  effects. 

Therefore,  with  the  Green’s  function  of  known  form  satisfying  all  boundary  conditions,  except  those 
on  the  failure  boundary  surface  itself,  we  can  obtain  an  integral  representation  of  the  dynamic  displace¬ 
ment,  u^d),  in  equation  (23.)  by  the  usual  methods.  That  is,  forming  the  inner  product  of  each  term  in  (23.) 
with  G™  and  integrating  over  the  source  coordinates  in  the  volume  V,  outside  the  failure  surface  and  over 
the  source  time,  gives: 


dV, 


—  V,  —  v,  —  V, 

Likewise,  forming  the  inner  product  of  each  term  in  equation  (28)  with  uk(d)  and  similarly  integrating  over 
the  source  variables  gives: 


Jdtjajp  wju^dV,  -  JdtJw  uk(d)dV,  =4nu«)(r  .  0 

—  V,  —  V, 


where  the  formal  integral  properties  of  the  della  functions  on  the  right  side  of  (28)  have  been  used.  Now, 
subtracting  the  second  of  these  equations  from  the  first  and  making  use  of  the  fact  that 


ultimately  gives: 


Tikd)d|Gkm  -  G  lk  3iUk<d) 
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The  first  integral  term  can  be  transformed  so  that  the  time  integration  can,  in  part,  be  evaluated. 
However,  it  is  clear  that  the  failure  zone  is  expanding  with  time  and  so  the  volume  V,  exterior  to  this  zone 
is  time  dependent.  Thus  we  can’t  simply  commute  the  time  derivative  with  the  volume  integration  in  the 
first  integral  term,  but  must  take  account  of  the  time  variation  of  the  spatial  integration.  To  do  so  we  can 
make  use  of  the  transport  theorem  for  the  case  of  a  boundary  that  moves  at  a  rate  that  is  not  equal  to  the 
particle  velocity  in  the  medium. 

This  case  has  been  treated  by  several  authors,  for  example  by  Eringcn  (1975)  and  by  Minster 
(1974),  the  latter  giving  several  detailed  derivations  of  the  result.  (Archambcau  and  Minster,  1978,  use 
the  theorem  in  a  similar  application  as  well.)  Specifically,  for  any  function  F  that  is  dependent  on  the 
deformation  within  the  medium  (ie.  dependent  on  the  Eulcrian  coordinates)  we  have: 


^JFdV-lf  dV't  I  FU"da' 

v,(tl)  v,(t,)  1  av,(tl) 


(31) 


where  U  is  the  velocity  vector  of  the  boundary  surface  and  n  is  the  normal  to  that  surface. 

Now  if  we  take  the  integrand  in  the  first  integral  in  (30)  to  be  the  function  F,  then  we  can  apply  (31 ) 
directly.  We  therefore  have,  for  this  first  integral  term: 


However  the  first  integral  on  the  right  side  of  this  equation  vanishes  identically,  that  is: 
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This  follows  from  the  fact  that  the  medium  is  at  rest  initially  so: 


and  since  Gkm  is  causal,  so  that 


au<d) 


=  0,  as 


h 


oo 


Gkm(r,t:  r.t,)  |„M  = 


aGkm(r,  t;  r,  t,) 


l,=-=o 


(Actually  the  above  is  true  for  all  t<t,,  but,  in  particular,  when  t,  becomes  infinite  then  it  is  true  for  all 
receiver  times  t.) 


Therefore,  using  these  results  in  (30)  and  further  applying  the  ordinary  form  of  Gauss’  Theorem  to 
the  second  integral  in  (30),  gives: 


This  result  is  one  form  of  the  integral  Green’s  function  representation  of  the  dynamic  radiation  field,  u£°, 
produced  by  a  growing  failure  zone  in  an  initially  stressed  medium.  The  result  is  formally  exact. 

The  first  integral  term  on  the  right  side  of  (32)  involves  an  integration  over  the  failure  boundary  and 
the  integrand  in  brackets  contains  the  same  combination  of  terms  as  are  involved  in  the  conservation  rela¬ 
tion  (25).  Thus,  the  integral  representation  contains  the  natural  boundary  conditions  that  apply  to  the 
physical  process.  Since  only  dynamical  field  quantities  are  involved  in  this  surface  integral  term,  it  is 
appropriate  to  interpret  this  integral  as  representing  the  interaction  of  the  dynamic  field  with  the  (grow¬ 
ing)  failure  boundary.  As  such  it  would  represent  scattering  and,  possibly,  absorption  effects  at  the  boun¬ 
dary.  (Note  that  the  rupture  rate,  UR  =  U,n|,  appears  in  the  integrand  of  the  integral  and  that  UR  is  a  func- 
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tion  of  the  energy,  L,  absorbed  in  the  process  of  failure.) 

The  second  volume  integral  contains  the  inertial  term  arising  from  changes  in  the  equilibrium  field 
in  the  elastic  medium  surrounding  the  failure  zone  and  corresponds  to  the  fundamental  source  term  that 
gives  rise  to  the  dynamical  response  of  the  medium.  Because  this  term  generally  represents  a  relaxation 
of  stress  and  a  reduction  of  stored  strain  energy  (which  shows  up  as  energy  in  the  radiated  elastic  wave 
field)  the  elastic  waves  produced  by  this  source  term  have  been  called  relaxation  source  fields  (Archam- 
bcau,  1964,  1968)  and  used  to  approximate  earthquake  radiation  fields.  This  term  represents  a  release  of 
energy  from  the  entire  prestressed  medium  surrounding  the  failure  zone  and,  as  such,  is  the  fundamental 
source  of  the  clastic  wave  radiation  produced  by  failure.  In  this  regard  the  gravity  term  appearing  with  the 
relaxation  source  term  involves  the  dynamic  changes  in  the  gravity  field  due  to  density  changes  in  the 
medium.  As  noted  earlier  and  as  described  in  the  Appendix  1,  it  is  very  small  compared  to  the  relaxation 
term  and  can  be  neglected. 


Vll.  Equivalences:  Initial  Value  (Relaxatiui.)  Sources 

As  was  just  mentioned,  the  result  given  in  (32)  is  slightly  different  than  that  obtained  in  earlier  work 
using  a  different  method  of  derivation,  which  was  based  on  the  view  that  the  radiation  process  associated 
with  failure  in  a  stressed  medium  could  be  described  as  an  initial  value  problem.  Nevertheless  the 
representation  in  (32)  can  be  transformed  to  a  form  that  is  nearly  identical  to  that  obtained  from  the  ear¬ 
lier  analysis. 

In  order  to  show  this  latter  "equivalence",  consider  the  integral  term  involving  the  inertial  source 
factor  (  3(i(p  3,  uk(e)’ ).  We  note  that  the  integrand  in  this  integral  may  be  "expanded"  to  the  form: 


at,  p  at ,< 


(e)' 


-P^eV(W) 


and  that  the  right  side  of  this  identity  can  be  used  in  the  integral  instead  of  the  quantity  on  the  left.  Thus 
the  integral  term  in  question  can  be  written  as: 


pa^Xcruv, 


dV, 


4! 


The  first  integral  term  on  the  right  side  of  this  equation  is  that  obtained  as  the  "relaxation  source"  term  in 
the  initial  value  formulation,  while  the  second  integral  is  an  extra  term  not  explicitly  present  in  the  initial 
value  formulation  results.  This  latter  term  can,  however,  be  recast  into  the  form  that  can  be  combined 


with  the  surface  integral  terms  in  (32)  by  using  the  transport  theorem  result  in  (31).  That  is,  using 
p  3t|uk(e)'G”  for  F  in  (31)  gives: 

4-00  -foo 

Jdt,  J  3ll[pdttu1[e)'(Gkro)]dV1  =  Jdt, 

—  V,(t.)  — 

(34) 


J  p  auuk(e)'GkmdV,  -  Jdt,  J  [p3liuk(e)'Gkm]u1n1dfl 

1  v.(i.)  —  av,(i1)L 


The  time  integration  in  the  first  integral  on  the  right  yields  the  spatial  volume  integral  evaluated  at  tlte 
(infinite)  time  limits.  Therefore  the  first  integral  on  the  right  is: 


j*. 

-t  J  pa^ofw, 

. 

J  pa^Gfdv, 

— *o 

v.d.) 

v.d.) 

(35) 


where  the  integral  vanishes  since  the  equilibrium  field  is  taken  to  change  continuously  (the  failure  zone 
growth  described  by  the  rupture  rate  UR  is  taken  to  be  a  continuous  function  of  time)  and  since: 

lim  9tluke)  =0  ;  limG^r,  t;  r,,  l,)  =  0 


Thus,  using  (35)  and  (34)  in  the  equation  (33)  gives: 


-Jdt,  J  3w[p  dfcU^GfdVj  =  Jdt,  J  p  d,,uk(e>  d,  G^dV ,  +  Jdt,  J  [p  3liU <e>Gkm] U.n.dfl , 


V,(i.) 


av,(t.) 


(36) 


Inserting  this  result  in  the  basic  representation  integral  for  the  dynamic  field,  equation  (32),  there¬ 
fore  produces: 


4jtu'm(r ,  l)  =  4rcu£)  +  Jdt,  J  { 


3iA 


^,  +  PVU 


GJ"  - 


dGkm 

Gff  +  P— U, 


uk(d)  } 


njd^7 , 


4-00 


4-00 


Jd''  I f 

—  V,(1.) 


8uk(c)' 


3t, 


3Gkm 


at, 


dv,  + 


Jdt,  J  pgk(d)Gk" 

—  v,(t.) 


dV, 


(37) 


where  a  factor  has  been  added  to  both  sides  of  the  result  in  order  to  put  the  final  equation  in  a  form 
comparable  to  that  given  by  Arehambcau  and  Minster  (1978).  In  this  result  we  have  defined  relative  dis¬ 
placements,  measured  from  the  initial  state,  so  that  as  previously  defined 


u(e)'  =  u(e)  -  ua) 


for  the  equilibrium  fields,  and  similarly 


u  =  u  -  u 


(I) 


for  the  "total"  fields  (dynamic  plus  equilibrium  fields).  Thus  relative  fields  are  measured  from  the  fixed 
initial  equilibrium  stale  (with  displacement  u0))  and  denoted  by  primes  throughout.  In  view  of  these 
definitions  and  the  relation  u  =  u(d)  +  u(c),  then 


u'  =  u(d)  +  u(c)' 


If  the  start  and  completion  of  the  failure  zone  is  taken  to  be  the  time  interval  (0,  x),  then  the  interval 
of  time  integrations  appearing  in  (37)  can  be  reduced.  That  is,  U,n,  =  0  and  3l|uk(e)  =  0  for  t,<0  and  for  t,>x 
in  the  integrals,  and  therefore  the  lime  limits  on  integrals  involving  these  factors  may  be  reduced  accord¬ 
ingly.  In  particular,  they  become  the  interval  (0,  min  [t,  x]),  as  noted  by  Archambeau  and  Minster  (1978). 
With  the  convention  that  the  failure  zone  development  occurs  over  the  time  interval  (0,  x),  then  we  also 
have  that  u(e)  =  ua)  -  ua)  =  0,  for  <  0;  while  u(e)  =  u(F)  -  u®,  for  t,  >  x,  with  u(F)  denoting  the  final  equili¬ 
brium  field. 

With  the  conventions  of  a  specific  reference  state  and  a  finite  source  time  interval  from  0  to  x  (which 
involve  no  loss  of  generality),  we  get  similar  results  from  (37)  as  those  given  by  Archambeau  and  Min¬ 
ster.  In  particular,  the  identical  volume  integrals  corresponding  to  relaxation  effects  and  gravitational 
effects.  However,  the  surface  integral  in  (37),  while  of  the  same  form  as  that  obtained  from  the  initial 


value  formulation,  is  nevertheless  different  in  detail  than  that  obtained  earlier.  In  particular,  the  dynamic 
stresses  and  displacements,  x£d)  34- x£e)  and  m  uk  -  uk(e),  appear  in  the  surface  integral  in  (37)  while  in 
the  initial  value  formulation  these  are  (implicitly)  replaced  by  the  "relative"  stress  and  displacement  com¬ 
ponents  tV  s  -  Xj$P  and  u'k  =  uk  -  uka).  The  fields  x^  and  u£e)  are  parametrically  time  varying  during  the 
failure  process  (by  definition)  and  the  fields  x£}  and  uk0>  are  fixed,  time  independent,  values  of  equilibrium 
stress  and  displacement  before  initiation  of  the  failure  process. 

Considering  the  simplicity  and  rigor  of  the  derivations  leading  to  the  representations  in  (32)  and 
subsequently  to  (37),  it  is  concluded  that  (37)  is  a  proper  exact  representation  of  an  initial  value  formula¬ 
tion  of  a  spontaneous  failure  source  (earthquake).  Thus  it  appears  that  the  direct  initial  value  foi.. relation 
obtained  earlier  did  not  express  the  surface  integral  term  completely,  or  at  least  that  this  latter  formulation 
was  imprecise  as  to  the  definitions  of  the  field  variables  in  the  integral  terms. 

Since  the  surface  integral  term  represents  scattering  and  absorption  of  energy  at  the  failure  zone 
boundary,  it  may  often  produce  smaller  second  order  effects  in  the  radiation  field  compared  to  the  relaxa¬ 
tion  integral  term  in  (32)  or  (37).  Indeed,  the  early  models  generated  from  the  initial  value  formulation 
entirely  neglected  the  surface  integral  term,  along  with  the  very  small  gravity  effects,  so  that  only  the 
relaxation  term  in  (37)  was  used  to  represent  the  radiation.  (The  resulting  field  predictions  were  called 
"transparent  source  approximations".)  Thus,  the  earlier  source  models  never  included  the  surface  integral 
contributions,  so  that  in  the  use  of  these  approximate  models  the  exact  nature  of  contributions  from  the 
surface  integral  term  never  arose..  (However,  neglect  of  the  surface  integral  contribution  was  done 
without  analytical  justification,  with  only  intuitive  arguments  used).  The  representations  in  (32)  or  (37) 
can,  in  any  case,  be  used  to  accurately  determine  the  nature  and  size  of  the  true  scattering-absorption 
effects. 

Some  concrete  justification  for  neglect  of  the  "scattering-absorption"  term  has,  however,  been  pro¬ 
vided  by  Stevens  (1981)  who  compared  the  "transparent  source  approximation"  fora  spherical  transition 
zone  in  a  prcstressed  medium  with  the  exact  solution  for  the  same  problem;  with  the  latter  containing  the 
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proper  surface  integral  scattering  term.  This  comparison  showed  that  the  surface  integral  did  indeed  con¬ 
tribute  small  reflection-refraction  effects,  but  that  a  much  larger  direct  first  arrival  radiation  field  was 
obtained  from  the  relaxation  integral  term  and  that  it  was  a  good  approximation  to  the  total  source  radia¬ 
tion  field.  However,  since  the  source  considered  by  Stevens  was  an  instantaneously  created  spherical 
failure  zone,  it  did  not  contain  the  possibilities  for  absorption  of  energy  as  would  a  spontaneous  failure 
process.  Thus  the  "absorption"  effects  at  a  growing  boundary  were  not  included  in  the  comparison. 
While  these  latter  effects  might  in  fact  "damp  out"  scattering  from  the  source  boundary,  so  that  the  two 
effects  might  cancel  each  other  out  in  part,  it  cannot  be  confidently  concluded  that  this  would  always  be 
the  case.  That  is,  it  is  not  clear  that  absorption  is  of  the  same  order  as  scattering  in  its  effect  on  the  total 
radiation  field  and  could,  at  least  in  some  circumstances  and  at  some  points  on  the  failure  surface,  be  a 
considerably  larger  effect  and  so  significantly  modify  the  total  field  from  that  predicted  by  the  relaxation 
term  alone.  (For  example  this  might  well  be  the  case  for  radiation  contributions  from  near  the  front  of  an 
advancing  thin  cllipsoidally  shaped  failure  zone). 

Therefore,  given  all  these  possibilities,  retention  of  the  surface  integral  contributions  in  specific 
source  models  seems  desirable  in  most  cases  of  spontaneous  failure.  Ideally  then,  models  of  earthquakes 
should  include  the  scattering-absorption  surface  integral  as  well  as  the  relaxation  term,  whether  (32)  or 
(37)  is  used  as  a  basis  for  the  model  predictions. 

It  is  also  appropriate  to  note  that  the  earlier  initial  value  development  employed  the  device  of  con¬ 
sidering  the  rupture  process  as  being  a  series  of  elemental  discontinuous  changes  in  the  rupture  dimen¬ 
sions,  with  these  changes  producing  corresponding  discontinuous  changes  in  the  equilibrium  field  uk(c>. 
By  taking  these  changes  to  be  of  infinitesimal  size,  and  the  time  interval  between  them  to  be  of 
infinitesimal  duration,  a  limiting  process  of  summing  the  infinitesimal  contributions  was  used  to  produce 
the  effects  of  continuous  failure  growth.  By  contrast,  the  present  approach  treats  rupture  growth  and 
changes  in  the  equilibrium  field  to  be  continuous  from  the  beginning.  These  two  methods  should  give  the 
same  results,  provided  the  failure  process  proceeds  in  a  continuous  fashion.  If,  however,  there  is  a  finite 
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discontinuous  change  in  the  rupture  rate  and  a  corresponding  discontinuity  in  the  equilibrium  field,  then, 
as  noted  by  Stevens  (1980,  1981),  the  integral  in  equation  (35)  will  not  be  zero.  However,  for  earth¬ 
quakes  which  result  from  spontaneous  rapid  failure  resulting  from  slow  loading,  it  is  not  likely  that 
discontinuous  changes  in  the  rupture  rate  occur,  so  that  the  imposition  of  a  continuous  time  varying  rup¬ 
ture  rate  in  the  formulation  seems  justified.  Nevertheless,  discontinuous  changes  cannot  be  ruled  out  with 
absolute  certainty.  Such  cases  can,  however,  be  considered  in  either  formulation  provided  care  is  taken  to 
explicitly  account  for  any  finite  and  instantaneous  changes  in  the  rupture  rate  and  the  equilibrium  field 

M(e) 

Uk  . 

An  example  of  a  singular  case  in  which  a  "supersonic"  rupture  rate  combined  with  spherical  sym¬ 
metry  of  the  developing  failure  zone  produces  discontinuous  behavior  occurs  when  an  explosion  is 
detonated  in  a  prestressed  medium.  In  this  case  the  failure  zone  can  be  completely  formed  by  a  super¬ 
sonic  shock  wave  before  the  medium  outside  this  failure  zone  can  react  dynamically,  since  the  shock 
velocity  (and  rupture  rate)  is  taken  to  be  larger,  at  all  points  on  the  spherical  failure  surface,  than  the 
highest  intrinsic  velocity  of  signal  propagation  in  the  medium.  Consequently  the  failure  zone  is,  in 
essence,  formed  instantaneously  in-so-far  as  the  (causal)  relaxation  of  the  equilibrium  field  external  to 
this  failure  zone  in  concerned.  This  type  of  source  has  been  treated  in  detail  by  Archambcau  (1968, 1972) 
as  an  initial  value  problem,  in  the  same  spirit  as  in  the  case  of  spontaneous  failure,  but  where  the  discon¬ 
tinuous  change  in  the  equilibrium  field  at  the  conclusion  of  the  shock  induced  failure  process  is  accounted 
for  explicitly  from  the  beginning.  Stevens  (1980,  1981)  obtains  the  same  result  using  u  method  similar  to 
the  "time  dependent  equilibrium  method"  described  here,  but  tailored  to  the  case  of  the  instantaneous 
spherical  failure  zone. 

In  the  present  formulation,  if  we  impose  the  conditions  of  this  singular  problem  on  the  integral 
representation  in  either  of  the  equivalent  forms  given  by  (32)  and  (37)  we  get  the  same  results  obtained  by 
Archambcau  and  Stevens;  even  though  the  representations  in  (32)  and  (37)  were  derived  under  the 
assumptions  of  continuous  failure  zone  growth  with  a  finite  rupture  rate.  In  particular,  taking  the  time  at 
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which  the  source  is  initiated  (ie.  t,  =0)  as  being  the  time  at  which  the  rupture  zone  has  been  formed,  so 
that  =  UR  =  0  for  all  times  equal  or  greater  than  t,  =  0,  then: 


u(e)'(r,  t,)H 


0  ;  t,  <  0 

ju(F)  -  u®]  H(t,) ;  t1  >  0 


where  H  (t,)  is  a  step  function  in  time  centered  at  t,  =  0.  Now  both  (32)  and  (37)  give: 


4ru^d)(r,  t)  =  J dt,  J  [tfGkm  -  G£  uk(d)] n,d Q ,  +  Jpu* 


o  av, 


i)Gk 


at, 


dV, 


(38) 


where  uk  = 


u^  -  uka) 


and  where  V!  is  the  fixed  volume  exterior  to  the  failure  zone.  This  is  the  same 


result  obtained  when  the  discontinuous  effects  arc  explicitly  accounted  for  when  deriving  the  integral 
representation.  Thus  the  results  in  (32)  and  (37)  apply  to  this  discontinuous  case  as  well  as  to  all  continu¬ 
ous  rupture  eases. 


VIII.  Equivalences:  "Dislocation-Stress  Pulse"  Sources 

Another  form  of  equivalent  representation  of  the  radiation  field  due  to  failure  in  a  stressed  medium 
can  be  derived  from  (32)  and  this  representation  involves  only  integrals  over  the  failure  boundary.  There¬ 
fore  it  can  be  compared  to  kincmatical  representations  (dislocation  models)  and  to  the  so-called  stress 
pulse  representations  that  involve  only  integrals  over  the  failure  surface.  The  approach  followed  is  a  sim¬ 
ple  generalization  of  one  used  by  Stevens  (1981),  which  was  introduced  to  demonstrate  a  similar 
equivalence  when  a  failure  zone  is  "instantaneously"  created.  In  this  case  we  want  a  result  similar  to 
Stevens’  result,  but  for  the  general  case  of  a  finite  rupture  rate. 

To  obtain  the  desired  representation  we  need  to  transform  the  relaxation  integral  term  in  (32)  to  a 
surface  integral  taken  over  the  failure  boundary.  To  do  so  we  will  need  to  develop  a  purely  formal 
integral  identity  involving  the  equilibrium  field  uk<c)'  and  the  dynamical  Greens  function  Gkm.  This  is 
achieved  by  considering  the  basic  differential  equations  satisfied  by  these  fields,  which  are: 
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-pgke)' 

d\G£  -  du(p  d,Gkm)  =  -  47tAkm(r -  r,;  t  -  t,) 


Taking  the  inner  product  of  the  equation  for  Gkm  with  on  both  sides,  and  integrating  over  the  source 
time  and  spatial  coordinates,  gives: 


+••  4*0  +®o 

Jdt,  J  [ati(p  d.Gnju^'dV,- Jdl,  J  3,G £  uk(c) dV ,  =  4rt  Jdt,  J  A^'dV,  =  4*u«’ 

—  v,(l.)  —  V,(t,)  —  V,(t,) 


Now  taking  the  inner  product  of  the  first  equation  with  Gkm  on  both  sides,  integrating  and  subtracting  the 
result  from  the  previous  equation,  gives: 


-H* 


+  jdt.  J  [pg^']GkmdV1=4J<>' 

—  V,(t,) 


(39) 


Noting  that 


G£a,uf>'=4*>'3,Gk"\ 


due  to  the  symmetry  of  the  clastic  tensor  Cljk|,  then  it  follows  that: 


a,G  £  uk(e)'  -  3,4e)'Gkm = 3,  [G  u  <e>'  -  x£e)'  Gkm 


Thus  the  integrand  in  the  second  integral  on  the  left  in  (39)  can  be  replaced  with  the  divergence  tenm  in 
the  equation  above  and  in  this  form  the  resulting  volume  integral  can  be  transformed  to  a  surface  integral 
using  Gauss’  theorem.  We  therefore  obtain  from  (39): 


Jdt,  J  [a,,(p  a,,Gkm)]  u^'dV ,  =  4jtu +  J dt,  J  ^G £  uke)  -  t£e) G™] n,dfl 

—  v,«,)L  —  av,o,)L 


-  Jdh  J  [pgke)' 

—  v,(i,) 


Gkmd  V , 


(40) 
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However,  we  note  that 


3jpa,Gr)]u<e)- 


d,,(p  at  uk(e>') 


Gkm  =  r\ 


pKg 


m.t  (c*V 


U 


-  p  a,  ,uk(c)'G 


m 

k 


Therefore  we  can  replace  the  volume  integral  term  on  the  left  side  of  (40)  using  this  last  equation.  This 
produces: 


Jdt,  J  [3li[Paiiukc>']]GkmdV1=4iruir,'+  Jdt,  J  [cj  uk">'  -  t^G^dV, 

—  V,{i,)  — «=  v,(*,) 


M 

—  V,(T.) 


_3_ 

at, 


(p  3,  Gkm)  uk(c)  -  (p  3.  uk(e) )  Gkm  dV 


M 


Pg^G^dV, 


v,(l.) 


In  addition,  we  can  apply  the  transport  theorem  to  the  volume  integral  involving  the  time  derivative  on 
the  right  side  of  this  equation  to  transform  it  to  a  surface  integral.  That  is,  using  the  transport  theorem 
given  in  equation  (31),  we  get: 


(p  ,Gk 


‘)uk(c)'  -  (p  3,uk<c))Gk" 


dV 


=\  J  [(p^GDu^'-fpa.u^r 


dV, 


V,(t,) 


Jdt,  J  [(P3liGkn,]uk(c)'-  [p3,,u<c)jGkm 


U,n,d£7 , 


av,(i,> 


where  the  first  term  on  the  right  vanishes  at  the  infinite  limits;  since  3, G™  =  Gkm  =  0,  for  t,-M-<»;  and 
uk(e)<  =  3tiuk(<:)’  =  0,  for  t,  — >  —  oo.  Therefore  the  previous  equation  can  be  rewritten  in  the  final  desired  form  as: 


Jdt,  J  [at,  [P  3liUk<e)] 

—  V,(i,) 


Gkmd  V ,  =4ku£y  + 


Jd‘,  1  j[c* 


+  PU,3lGkm 


av,d,)t 


uk 


(c)' 


hi0'  +  P  G,  3,,uk 


(e)' 


GkmMfl, 


Jdt,  J  pg^'c; 

—  ViUl) 


dV, 


(41) 


This  integral  identity  now  expresses  the  relaxation  term  that  appears  in  (32)  in  terms  of  a  surface  integral 
on  the  failure  boundary.  Thus  (41)  can  be  used  to  transfomi  the  basic  representation  integral  equation  in 
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(32)  into  the  following  equivalent  representation: 


4nu'm(r,  t)=  Jdt,  I  | It*  +  (p 3, u'^Ujlcr -  [c?  +  (p 3, .G^U, 
—  av,d,)lL 


+  Jdt,  J  p  g'kGk" 

—  v,(tl) 


dV, 


(42) 


where  we  have  a  result  expressed  in  terms  of  the  "relative"  displacement,  stress  and  gravity  fields;  that  is 
uk  =  uk(d)  +  uk(e),  etc. 

Equation  (42)  has  a  fomi  that  is  quite  similar  to  the  standard  elastodynamic  integral  equation  used 
to  obtain  equivalent  dislocation  "models"  for  earthquakes.  (If  U,n,  =  UR  were  set  to  zero  in  (42),  then  it 
would  have  the  identical  form.)  However,  the  integral  representation  in  (42)  is  actually  considerably  more 
complex  than  the  standard  form;  in  spite  of  its  relatively  simple  expression. 

In  order  to  display  this  complexity  and  isolate  the  various  contributions  to  the  radiation  field,  it  is 
better  to  rewrite  (42)  with  the  dynamic  and  equilibrium  field  explicitly  displayed.  That  is,  (42)  can  be 
rewritten  in  expanded  form  as: 


4rcu'm(r  ,  t)  =  J”  dt,  J  •  [t^  +  (p  3ti  ukd))U,j  Ckm  -  [g  J*  +  (p  3,,  Gkm)U,]  uk(d)j  n,  d£7 , 


av.ti.) 


Jdt,  J  |  +  (p  3ti  uke)  )U,j  G™  -  |g  *  +  (p  3ti  G  ™)U,j  uk<e)l  n,dfl ,  +  Jdt,  J  (pg'k)CkmdV,  (43) 

—  av,(i,>  J  v,do 


where  the  (time  varying)  equilibrium  field  contributions  are  explicitly  displayed  along  with  the  purely 
dynamical  parts  of  the  displacement-stress  fields.  Here,  of  course,  u'm  =  u^d)  +  u£r  on  the  left  hand  side  of 
the  equation,  so  that  the  relative  displacement  field  is  represented. 

The  boundary  conditions  at  the  failure  surface  are  such  that  the  rupture  rate  vector  U  is  a  function  of 
both  the  equilibrium  and  dynamic  fields  (see  equation  26).  Further,  from  equations  (18)  and  (25),  we 
have  the  boundary  conditions: 


fljuf’  nj  =  0  ;  re  3V, 
pU.ffvJ  =  -  M  ;  re 3V, 


(44) 


with  UR  =  U,nlt  corresponding  to  the  rupture  rate.  Here  vk  and  ^  =  r^n,  arc  the  total  velocity  and  traction 
vector  components,  involving  the  sum  of  equilibrium  and  dynamic  fields.  Consequently,  using  these  con¬ 
ditions  together  gives: 


pURH>  J  =  -  M 

Thus,  on  the  failure  boundary  3V,  (t,): 

PUrvV  +  i^pIV'V,  +  <V> 


(45a) 


(45b) 


where  °Vk  and  denote  field  variables  from  within  the  failure  zone.  Also,  the  velocity  fields  both 
inside  and  outside  the  failure  zone  have  the  form  v'k  =  v^d)  +  vke)‘ ’.  Thus  the  boundary  condition  applying 
on  the  failure  surface  involves  both  the  equilibrium  and  dynamic  velocity  components  in  the  particular 
combination  expressed  by  (45).  Therefore,  it  is  appropriate  to  rearrange  the  terms  in  (43)  so  that  the 
boundary  conditions  expressed  in  (45)  apply  to  one  of  the  surface  integrals,  while  the  first  of  the  condi¬ 
tions  in  (44)  applies  to  the  oilier  surface  integral.  That  is,  rearranging  (43)  to  accommodate  the  boundary 
conditions  in  separate  integral  terms  yields: 


Gkm  - 


gk”  uk  +  (PUR  3(i  G™)  u\ 


4?ru'm(r  ,  •)  =  J  dt,  J  |  ^(d)  +  pUR  3ti  u'k  ( 

—  av.d.)  I 

+oo  +"> 

+  J  dt,  J  [ lk<e)' Gkm  -  gkm  uke)j  d(2,  +  J  dt,  J  (p  g\)  Gkm  dV , 


dfl, 


(46) 


av,d,) 


where  gkm  denotes  the  tractions  associated  with  the  Greens  function  G".  In  this  form  the  second  surface 
integral  involves  only  equilibrium  field  changes  and  the  first  boundary  condition  in  (44),  pertaining  to  the 
equilibrium  field,  applies  to  it.  In  this  regard  we  note  that  the  first  boundary  condition  in  (44)  applies  to 

the  equilibrium  traction  t/kc)  n,  =  (ifc ’ '  +  x^)n,.  But  since  [[  n,  ]]()  also,  then  it  follows  that  traction 
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continuity  applies  to  the  relative  stresses  as  well.  That  is 


Itf  nJ  =  0 

and  this  relation  applies  in  the  second  integral.  On  the  other  hand,  the  first  surface  integral  involves 
dynamic  and  equilibrium  field  changes  at  the  boundary,  with  the  boundary  condition  (45)  applying. 

If  there  were  no  equilibrium  field  changes,  then  the  second  integral  in  (46)  would  vanish,  since  the 
equilibrium  stress  and  displacement  fields  appearing  in  (46)  are  measured  relative  to  the  initial  state  of  the 
medium.  In  this  case  there  could  be  no  failure  boundary  growth,  so  UR  =  0  is  implied,  and  the  first 
integral  would  reduce  to  the  standard  form  with  V!  and  3V,  fixed.  Then  we  would  have: 

-foo  +*»» 

47uiT  =  J  dt,  J  Gkm  -  gkm  u,<d>]  d(2 ,  +  J  dt,  J  p  gk<d)  Gkm  dV, 

—  av, L  —  v, 


and  the  first  integral  would  represent  pure  scattering  from  the  fixed  surface  3V,,  while  the  volume  integral 
would  be  non-zero  if  there  were  dynamical  changes  in  the  gravity  field. 

Using  this  limiting  case  as  a  guide  to  the  interpretation  of  (46),  we  can  rewrite  this  equation  as: 


47tU 


n(r,t)=Jdt,  |  [t<d)Gkm-gkmuf]da,  +  {dt,  J  pUR[(3liu'k)Gkm-(31,Gkni)u'k]da, 

+«*>  -foo 

+  J  dt,  J  [t^  Gkm  -  gkm  u-']  Aa ,  +  J  dt,  J  p  g'k  Gkm  dV, 


(47) 


av,(t,) 


and  interpret  the  first  integral  as  the  term  due  to  scattering,  the  second  integral  as  that  due  to  absorption  of 
energy  along  the  growing  boundary  (because  of  the  presence  of  UR  and  its  dependence  on  the  energy,  L. 
required  for  failure)  and  the  third  integral  as  being  the  equivalent  for  stress  relaxation  and  energy  release 
in  response  to  the  growing  failure  boundary.  Here  again,  the  last  integral  is  the  (generally  negligible) 
gravity  term.  Because  of  the  boundary  condition  in  (45),  the  first  and  second  integrais  arc  actually  con¬ 
nected  -  which  would  be  expected  since  scattering  and  absorption  at  the  failure  boundary  must  be  inter¬ 
connected,  as  was  noted  earlier. 
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The  field  decomposition  can  be  symbolically  expressed  as: 


4n  u'  (r  ,  t)  = 


.  (u(d))  +  I,  (u(d)  +  uw') 


+  IR  (uw  )  +  L(g') 


(48) 


where  the  terms  on  the  right  denote  the  scattering  (Isc),  absorption  (I,),  relaxation  (IR)  and  gravity  (Ig) 
terms  in  (47).  Here  only  the  dependence  on  the  dynamic  and/or  equilibrium  displacements  arc  explicitly 
displayed,  since  the  corresponding  tractions  appear  with  these  displacements  in  the  integral  terms  and  can 
be  derived  from  them.  (The  brackets  arc  used  to  indicate  that  the  two  integral  terms  I5C  and  la  should  be 
treated  together  in  solving  this  equation.) 

The  relaxation  term  in  equation  (47)  is  now  in  the  form  of  a  surface  integral,  with  the  possibility 
that  it  can  be  viewed  in  terms  of  equivalent  stress-pulse  and  dislocation  sources.  In  particular,  we  can 
write  this  term  as: 


Ir  (u(e) )  =  4n  u®  +  u?» 


(49a) 


with  u®  and  as  displacement  fields  due  to  equivalent  stress-pulse  and  dislocation  type  sources, 
defined  by 


4mi'S)=Jdt1  J  t^Gfdfl, 

—  av.d.) 

+oo 

4mii?)s-Jdti  J  u^gTdfl, 

—  av,(tl) 


(49b) 


These  two  equivalents  are  not  independent  however,  since  once  the  equilibrium  displacement  change  uk(cV 
has  been  determined  for  the  static  inclusion  problem  then  tk(c)’  may  be  simply  derived  from  it.  Further,  for 
the  static  inclusion  problem  involved  here,  we  have 


uk(e)-uk(e) 


+  uk' 


(i) 


where  ukc)  and  u™  arc  the  displacement  fields  with  and  without  a  failure  zone,  while  uk(c)  represents  the 
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change  in  the  initial  state  due  to  the  presence  of  the  failure  zone.  Here,  uk(c)  u^'  and  uk0)  arc  all  com¬ 
ponents  of  biharmonic  vector  fields.  Thus,  we  already  know  the  spatial  form  of  both  uk(c!  and  tk(c\  and  that 

/  / 

the  coefficients  of  the  biharmonic  series  for  uk(c) ,  and  the  corresponding  traction  t^e) ,  are  parametrically 
dependent  on  the  source  time  c,  due  to  the  variable  dimensions  of  the  failure  zone  with  time.  Therefore 
the  surface  integrals  in  (49)  can  be  analytically  evaluated  once  the  failure  surface  is  specified. 

These  equivalent  source  integrals  can  be  expressed  in  more  explicitly  recognizable  form  if  we  use 
the  definition  of  uk(e)'  in  terms  of  field  changes  from  the  initial  state.  That  is,  we  can  also  write  (49)  as: 


These  two  source  equivalents,  corresponding  in  sum  to  the  relaxation  source  field  term,  may  be  inter¬ 
preted  as  a  continuously  variable  "stress-pulse"  on  the  expanding  failure  surface,  with  magnitude  equal  to 
the  traction  "drop”  on  this  surface,  plus  a  continuously  expanding  closed  dislocation  surface  with  a  dislo¬ 
cation  "offset"  equal  to  the  displacement  change  from  the  initial  state  at  the  failure  boundary. 

Thus  we  find  that  failure  induced  volumetric  stress  relaxation  and  the  seismic  radiation  associated 
with  it  can  be  accounted  for  by  a  superposition  of  two  familiar  surface  distributed  source  equivalents; 
namely  a  pair  of  related  stress-pulse  and  dislocation  equivalents.  Since  the  radiation  fields  produced  by 
each  of  these  equivalents  will  be  superposed  on  one  another,  one  can  expect  a  total  field  having  charac¬ 
teristics  that  may  be  quite  different  from  cither  of  the  equivalents  separately.  In  any  case  the  equivalent 
form  of  the  integral  equation  (47),  or  alternately  as  expressed  in  (48)  and  (49),  constitutes  a  second 
representational  form  for  the  clastodynamic  radiation  produced  by  failure  of  a  stressed  medium  and  we 
sec  that  this  equivalent  form  contains  not  only  combined  stress-pulse  and  dislocation  equivalents  that 
account  for  volume  relaxation  effects,  but  boundary  absorption-scattering  terms  that  arc  not  usually 
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negligible.  Therefore  it  is  generally  necessary  to  solve  the  integral  equation  in  (47)-(48)  to  obtain  the 
total  wave  radiation  field.  These  applications  arc  investigated  in  detail  by  Archambcau  and  Dills  (1989). 

IX.  Summary  and  Conclusions 

Examination  of  some  of  the  previously  developed  representations  for  seismic  radiation  from  earth¬ 
quakes,  particularly  those  using  the  divergence  of  the  medium  stress,  or  the  "stress  glut",  in  the  moment 
tensor  representatioas  of  earthquake  radiation  fields  leads  to  the  conclusion  that  these  representations 
have  no  rational  physical  basis,  nor  can  they  be  logically  related,  in  even  an  approximate  sense,  to  either 
the  kinematics  or  dynamics  of  such  a  source. 

The  appropriate  representation  of  the  radiation  field  can,  however,  be  obtained  in  a  quite  straightfor¬ 
ward  manner  by  noting  that  the  displacement-stress  fields  in  the  medium  are  a  sum  of  dynamic  and  equili¬ 
brium  components  and  that  the  equilibrium  field  component  is  lime  dependent  when  there  is  spontaneous 
failure  within  the  medium.  It  is  shown  that  this  temporal  variation  in  the  equilibrium  field,  which  is  due 
to  creation  of  a  failure  boundary  within  an  initially  stressed  medium,  gives  rise  to  an  "equivalent  force" 
term  in  the  equations  of  motion  in  the  linear  zone  outside  the  failure  boundary  and  corresponds  to  the 
inertial  effect  of  changes  in  the  equilibrium  displacement  field  in  the  medium  surrounding  the  failure 
zone.  Further,  it  is  shown  that  this  effect  only  occurs  when  the  medium  is  prestressed  and  a  new  boun¬ 
dary,  enclosing  a  zone  of  material  altered  by  the  failure  process,  is  created.  This  inertial  term,  rather  than 
a  term  involving  the  divergence  of  the  stress  in  the  medium,  is  the  proper  "equivalent  force  term"  associ¬ 
ated  with  an  earthquake.  However,  because  the  process  of  failure  is  spontaneous  and  depends  on  the 
dynamic  radiation  field  for  its  continued  growth,  the  near  field  radiation  effects  arc  complicated  by  energy 
absoiption  and  failure  boundary  scattering  that  will  also  manifest  themselves  in  the  total  radiation  field 
observed.  These  latter  effects  are  determined  by  the  boundary  conditions,  expressing  conservation  laws, 
on  the  growing  failure  surface. 

In  order  to  treat  both  the  inertial  or  relaxation  effects  in  the  medium  surrounding  the  failure  boun¬ 
dary  and  the  boundary  scattering  and  energy  absorption  associated  with  failure  growth,  it  is  appropriate  to 
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reformulate  the  representation  of  the  problem  as  an  integral  equation  involving  both  the  relaxation  effects 
and  the  boundary  interaction  effects.  A  Green’s  function  integral  equation  for  the  displacement  in  the 
linear  region  outside  the  failure  zone  volume  is  therefore  developed  using  both  the  equations  of  motion 
and  the  boundary  conditions  on  the  (growing)  failure  surface.  The  resulting  integral  representation  of  the 
radiation  field  involves  a  volume  integral  accounting  for  the  relaxation  of  the  equilibrium  field  in  the  elas¬ 
tic  medium  surrounding  the  failure  zone  (the  dominant  effect)  and  surface  integral  terms  that  may  be  indi¬ 
vidually  identified  as  being  associated  with  scattering  and  energy  absorption.  Further  it  is  shown  that  this 
"representation  theorem"  reduces  to  the  classical  result  in  cases  when  there  is  no  failure  at  all,  or  when  the 
failure  process  has  stopped. 

The  integral  representation  obtained  is  shown  to  be  susceptible  to  transformation  to  other  forms, 
which  is  useful  from  a  computational  point  of  view  as  well  as  being  important  from  a  conceptual  stand¬ 
point.  In  particular,  it  is  shown  that  the  integral  term  representing  the  relaxation  related  radiation  can  be 
put  in  a  form  identical  to  that  obtained  in  previous  work  (eg.  Archambeau,  1968;  Archambeau  and  Min¬ 
ster,  1978)  which  used  an  initial  value  approach  to  generate  the  representation  of  the  radiation  field. 
However,  it  is  found  that  the  surface  integral  terms  representing  scattering  and  absorption  at  the  failure 
boundary  are  more  precisely  defined  by  the  representation  obtained  here  and  have  a  slightly  different 
form  that  represents,  at  the  least,  a  clarification  of  the  earlier  results.  Alternately  the  volume  integral 
representing  seismic  radiation  from  relaxation  in  the  medium  surrounding  the  failure  zone  can  be 
transformed  into  surface  integrals  over  the  failure  zone  boundary  which  can  be  identified  as  the  superposi¬ 
tion  (or  sum)  of  a  surface  distributed  dislocation  and  a  distributed  stress-pulse  involving  only  the  time 
dependent  changes  in  the  equilibrium  fields  at  the  boundary.  These  two  equivalents  are  not  independent 
however,  since  they  arc  both  the  result  of  changes  in  the  equilibrium  field  from  its  initial  state.  (This 
latter  field  can  be  obtained  by  standard  methods  of  solution  for  inclusion  problems).  The  remaining  sur¬ 
face  integral  terms,  involving  both  the  dynamic  and  equilibrium  field  changes,  can  again  be  individually 
identified  as  terms  representing  scattering  and  absorption  at  the  failure  boundary. 
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The  integral  equation  representations  describing  the  failure  induced  seismic  radiation,  whether  those 
containing  a  volume  relaxation  term  explicitly  or  that  representation  with  this  particular  term  represented 
by  the  dislocation -stress  pulse  equivalent,  arc  exact  and  applicable  to  any  prestress  state  and  source 
geometry.  Further,  the  rupture  geometry  and  growth  rate  arc  prescribed  in  the  integral  representation 
vhen  the  boundary  condition  involving  energy  conservation  is  used. 

In  some  cases  it  is  likely  that  the  surface  distributed  dislocation-stress  pulse  equivalent  form  will  be 
more  convenient  from  a  computational  standpoint,  but  in  any  case  the  relaxation  integral  or  the  equivalent 
dislocation-stress  pulse  integrals  can  be  evaluated  independently  from  the  remainder  of  the  integral  equa¬ 
tion  if  a  failure  zone  geometry  is  assumed,  since  these  integrals  depend  only  on  the  equilibrium  field 
changes  which  can  be  obtained  independently  through  solutions  of  static  inclusion  problems.  This  means 
that  a  moment  type  expansion  can  be  obtained  for  the  main  (relaxation  term)  contribution  to  the  radiation 
field  by  standard  methods  and  that  corrections,  or  perturbations,  to  these  moment  terms  can  be  obtained 
by  iterative  approximation  of  the  scattering-absorption  integrals.  (These  approximations  are  considered 
in  detail  in  a  companion  study  by  Archambcau  and  Dills,  1989). 

Thus  the  integral  representations  obtained  here  form  the  basis  for  a  moment-tenser  representation  of 
earthquake  source  that  is  dynamical  in  nature  and  satisfies  the  required  conservation  laws  at  the  failure 
boundary,  as  \vc!i  as  iu  uii.  ..uudirg  lincar  acne.  Comparing  these  results  with  previous  representa¬ 
tions  of  earthquake  sources,  in  particular  the  "stress  drop",  "stress  glut",  "boundary  stress  pulse”  or  "dislo¬ 
cation"  representations,  leads  us  to  conclude  that  none  of  these  "models"  is  generally  appropriate, 
although  the  latter  two  may  be  adequate  approximations  in  certain  circumstances  since  they  may,  for 
some  failure  modes,  approximate  the  appropriate  dislocation-stress  pulse  equivalent  defined  here. 
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Appendix  1  -  Magnitude  of  Seismic  Radiation  Effects  Due  to  Induced  Gravity  Field  Changes 


In  the  initial  state  prior  to  failure,  mechanical  equilibrium  requires: 


(A-l.) 


where  the  gravitational  potential  <t>!  (x)  and  initial  density  p,  (x)  satisfy 


V2  =  4  7i  G  P[ 


(A-2a.) 


or, 


=  4  jr  G 


I 


Pi 

Ix-x'l 


dV 


(A-2h.) 


Here  rotational  forces  for  the  planetary  body  are  neglected  as  small  relative  to  gravity. 

The  final  state  of  the  medium  will  be  different  after  the  earthquake  but  will  also  satisfy  equilibrium 
conditions  and  (changed)  boundary  conditions,  so  that 

=  Pf  a, 

V2  d>F  =  4  rr  G  pF  [  '' 


and  therefore 


ch00  -  p,  d,  d>0) 


Noting  that  the  left  side  of  this  equation  involves  the  stress  drop,  then 


(A-4.) 

where  tf  is  the  spatial  source  term  used  by  Gilbert  to  represent  the  seismic  source  "associated"  with  an 
earthquake. 


Now  setting 


so  that 


Pf=Pi  +  8p  1 
<I>f  =  <&i  + 


(A-5.) 


V2  ct)F  =  V2  <t>,  4  V2  (60>)  =  4  n  G  (p„  4  8p) 


59 


then  we  have  immediately: 


V2  (SO)  =  4  Ji  G  (8p) 


(A-6a.) 


Therefore  the  change  in  the  gravitational  potential  is  given  by 


80  =  4  G 


j-rKr  dv 

i  i*-*  i 


(A-6b.) 


Further,  by  the  equation  of  continuity 


5p  =  -  V  •  (p,  s) 


(A-7.) 


where  s  denotes  the  displacement  field  produced  by  the  earthquake. 

Inserting  (A-5.)  in  (A-4.),  and  using  (A-6b.)  and  (A-7.)  in  the  result,  gives: 


r  V  •  (p,  s)  , 

y'“-4"0  p'v[|-i^TdV 


+  V  •  (p,  s)  V 


f_£L_ 


dV  +  V  •  (pi  s) ' 


jrdV 


(A-8.) 


Since  we  arc  interested  in  the  order  of  magnitude  of  the  source  term  /,  then  the  initial  density  can  be 


approximated  as  a  constant,  p,  in  the  integrals  and  in  this  case  (A-8.)  has  the  approximate  value: 
-4  ]tGp1{(l  +  V-s)  V  jj  dvj  +  (V-s)  V  jj 


(A-9). 


Noting  that  the  divergence  of  the  displacement  field  is  the  dilatation  and  that  it  is  small  compared  to 
unity,  that  is  observationally  of  maximum  value  near  1CT3  for  earthquakes  at  near  distances,  then  the  fac¬ 
tor  (1  +  V  •  s)  in  the  first  term  in  (A-9.)  can  be  approximated  by  unity.  Therefore: 


7V=  -4  jiG  p2 


V dV  +<?-s)v 

I X  -  X  | 


dV 

Ix-x'l 


(A- 10.) 


To  this  point  the  results  are  general,  in  that  the  volume  of  integration  for  the  gravitational  effect  can 
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be  taken  as  the  entire  volume  exterior  to  the  failure  zone  (the  volume  denoted  V,  in  the  text)  or  the 
volume  interior  to  the  failure  zone  boundary  denoted  V'  in  the  text,  or  both.  Gilbert’s  representation 
involves  an  integration  over  V,,  while  that  of  Backus  and  Mulcahy  involves  an  integration  over  V'. 

The  representation  of  the  force  equivalent,  f,  used  by  Gilbert  is  given  by: 

f  =  yv  H(t)  ;  r'e  V, 

where  H(t)  is  a  step  function.  The  designation  of  the  region  of  applicability,  r'e  V,,  means  that  this 
equivalent  is  defined  over  the  volume  exterior  to  the  failure  zone.  On  the  other  hand,  the  "stress  glut"  for¬ 
mation  of  Backus  and  Mulcahy  has  this  same  form  in  the  low  frequency  limit,  but  with  the  region  of 
applicability  for  the  equivalent  defined  to  be  over  V',  the  failure  zone. 

Considering  the  case  appropriate  to  Gilberts  formulation,  for  which  (A-10.)  is  directly  applicable, 
and  noting  that 

Yto  (6  ,  <t>)  Yj™  (6',  O')  ;  r  >  r' 

Y|^  (0  ,  <t>)  Yta  (9',4>')  :  r  <  r' 

where  Ylm  (0  ,  c»  =  P|m(cos  0)  cim®,  with  P™  the  associated  Legendre  function  and  where  Yb*n,  is  the  complex 
conjugate  of  Y^  ,  then  in  (A-10.)  we  have: 


v 

f  ^ 

4k 

1 

l,m 

21+1 

*.  * 

. 

Ix-x'l 

£ 

l.m 

4k 

21+1 

f  dV 

~  Y7 

r  R 

^  f  r'2  jt>  ,  a—  f  r'  jr' 

J  Ix-x'l 

—  V 

i  ur  +  4K  |  r  or 

r  0  r 

-4k/ 3  re, 


where  er  is  the  unit  vector  in  the  radial  direction.  Here  the  integration  volume  V,  has  been  taken  to  extend 
from  the  origin  to  a  large  distance  R  horn  the  failure  zone,  as  required  by  the  finite  dimensions  of  the 
planet.  Similarly,  the  integration  has  been  taken  through  the  failure  zone  in  order  to  give  an  upper  esti¬ 
mate  for  the  size  of  this  term. 

Likewise  the  magnitude  of  the  other  integrals  term  in  (A-10.)  can  be  estimated  by  noting  that  the 
dilatation  0  =  V  •  s  is  a  hannonic  function  (eg.  Love,  1944).  For  the  region  outside  the  failure  zone  litis 
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harmonic  function  is  always  of  the  form 


where 


1 

r"*1 


s.<e.4>) 


(A-l  1.) 


sn(0  .<)>)=  2 


Onn,  COS  m  4)  +  sin  m  <}) 


)  Pnm(cos  0) 


with  the  constants  ,  p„m  dependent  on  the  failure  zone  geometry  and  the  initial  stress  state  of  the 
medium.  (See  Archambcau,  1964  and  1968  for  examples.)  Thus,  with  |x-x'|  =  r*,  the  second  integral 
term  in  (A- 10.)  can  be  expressed  as: 


l 


V' s  j3_» 


I  X  -  X  | 


' 

'SnO',4.')' 

Jr 

«. 

j.'n+l 

r'2  sin  0'  d0'  d<|)'  dr' 


Now,  employing  a  method  due  to  Love  (1944),  and  applied  in  similar  circumstances  by  Archambeau 
(1968,  p.  255)  the  integration  over  V,  can  be  rearranged  lo  be  performed  over  surfaces  of  constant  r* 
around  the  "observation  point"  at  r  and  then  (f.nally)  over  r*.  to  cover  all  of  V,.  We  therefore  obtain, 
after  extending  the  region  of  integration  to  include  the  failure  zone  region  as  well  as  the  surrounding 
medium: 


I 


S,(9',<t>') 

(r'T1 


ACT 


where  R  >  r  and  Q*  is  the  spherical  surface  of  radius  r*  =  |r  -  r'|.  Since 


I. 


Sn(9',40 

(rT+1 


d  a’=< 


4ji  r 


Sn(9  ,  <t>) 


rn+1 

0 ;  r*  >  r 


;  r  <  r 


by  an  extension  of  the  mean  value  theorem  for  harmonic  functions  (sec  Archambcau,  1964  for  details). 
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then 


r 

f  d3*' =  471 S 

A  lX“Xl  » 

j  ,•  i,- 

0 

r”+ 


Consequently  the  second  term  in  (A- 10.)  becomes: 


1 


V  •  s  j3  , 
- d  x 


|x-x'| 


=  2k  Vfr2  ©)  =  2k  r2  V0  +  4  jt  r  0  er 


Now,  collecting  results  from  the  integral  evaluations  one  has  for  the  equivalent  force  in  (A- 10.): 


Yv  =  -  8n  G  p2 


4/3  r  0  er  + 


(A- 12.) 


To  directly  estimate  the  contribution  of  such  a  source  term  to  the  wave  field  we  observe  that  the 
dynamic  displacement  field,  u,  is  given  by: 


Uk(x  .  0  =  -j-  [  dto  f  f|  (x„  ,  U  G,k  (x  ,  x0  ;  t„  -  t)  d3x0 
4n  Jo  { 

with  f,  =  y,v  (x0)  H(t„)  and  G*  a  tensor  Greens  function  for  the  medium.  Alternately,  the  dilatation  and  rota¬ 
tion  of  the  dynamic  field  is  given,  in  the  spectral  domain,  by  the  simpler  relations  (eg.  Archambcau, 
1968): 


Xo(x.<o)=— ~  f  H  (Oc  i<l>l'dt0  f  (V  •  f)  d3x0 

471  Pvp  i  i  r 

X(x.a»)  =  — l—  f  H(gc"10>udto  f  {(Vxf)^-d3x„ 
4jtpv,2J  l  2  r 


(A- 13.) 


with  infinite  space  scalar  Greens  functions  used,  where  x‘  =  1  r-r0  |  and  kT,  =  w/vp,k,  =  u)/v5.  Here 
denotes  die  dilatation  while  x  =  (Xi  .X2.X3)  denotes  the  rotation  field  with  Cartesian  components  xr 
These  physical  potentials  should  be  of  the  order  of  the  strain  changes  observed  from  an  earthquake  if  this 
source  term  is  the  proper  equivalent  for  an  cartnquake  source.  These  potentials  arc  therefore  convenient 
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for  estimation  purposes. 


From  (A-12.)  we  have: 


V  x  f  =  ~  n2  P2  G 


V  •  yv  =  -  32  7U  2  p2  G 
r  V  0  x  e,j  = 


,,  5  3  0  1  r 

re+6rir+3e 


16rt2 


p2G 


a©  . _ i_  ae  . 

39  e*  sin  0  3<p  69 


where  8  has  the  form  given  in  (A-l  1.).  Inserting  these  relations  along  with  the  expansion  for  9  into  (A- 
13.)  and  introducing  the  spherical  wave  the  expansions  for  the  scalar  Greens  functions,  of  the  form 


-ikr 


=  -  ik  £  (2  1  +  1)  P[  (cos  y)  ji  (k  r')  h,(2)  (k  r)  ;  r  >  r' 
r  1=0 


gives: 


X*G)(r,ai)  = 


327I2  pG 
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,5  1 , 
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jn(kpr')  dr' 


Z(G)  (r ,  o))  =  -Z&p-  ZSn(0 

n  q 


n-1 

y  j„(k,  r')  dr' 


Here  S n(0 , 0)  denotes  a  vector  with  components  of  the  same  type  as  Sn(0  ,  p).  In  deriving  these  exprcs- 
sioas  we  have  evaluated  the  transform  of  the  step  function  H(t)  and  have  used  the  orthogonality  relation 


2k  it 


f  f  P,  (cos  Y)  sn  (0' ,  40  sin  9'  d0'  dp'  =  S„  (0  ,  p)  5,n  , 

00  2n  +  1 


with  y  the  angle  betweend  the  vectors  r  and  r'.  Again  the  range  of  integraton  has  been  extended  through 
the  failure  zone  to  the  origin,  in  order  to  obtain  maximum  estimates. 

In  evaluating  the  radial  integrals  in  ihc  expressions  for  the  wave  field  induced  by  gravity  changes 
we  can  further  simplify  results  and  obtain  an  accurate  estimate  by  taking  the  upper  limit  (R)  to  be  very 
large,  or  infinite,  and  using  the  integral  relation: 


J  x^"1  Jv  (a  x)  dx  -  2^  1  a  4 
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r(i/2v+  1/2  p) 
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valid  for  -  Re  u  <  Re  p  <  3/2.  Wc  have  in  this  case;  for  n  >  2; 


-(G)/  .  32jt2  p  G  „ 

Xo  («•.“)<  — ~  E 


n  k. 


n-3 


2n  r(n) 
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2n  T(n  +  1/2) 


S„(0  ,  <D)  h‘  '  (kpr) 


(A-14a.) 
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(A-14b.) 


2T(n+  1/2) 

The  magnitude  of  the  dilatation  is  larger  than,  but  of  the  same  order  as,  the  components  of  the  rota¬ 
tion  vector  field.  Therefore  it  is  sufficient  to  consider  the  si/c  of  only  Xcf’\  >n  (A- 1 4a.),  relative  to 
observed  dynamic  strain  changhes  accompanying  earthquakes.  To  obtain  a  numerical  estimate  of  Xu'\  it 
is  sufficient  to  consider  the  (extreme)  case  in  which  a  fluidized  spherical  failure  volume  is  produced  in  an 
initially  shear  stressed  medium.  In  this  case  wc  have  (Archambeau,  1968;  Landau  and  Lifshitz,  1959)  the 
specific  form  for  ©  in  (A-l  1.); 

©(r)  =  ~  S2  (9  ,  <t>) 


with 


S2(0  ,  <t>)  =  E  a2w  cos  m<f>  +  sin  mp  j  P£\cos  0  ) 

m=0  ^ 


where 


(Q2m)  =  p(7-"5a)'  R°3(  °  °  }  =  -2 


(P2m)=  p(7'4o)  R<’(  °  t,<2,/2  }  :  m  =  0’  1  ’2 

Here  x™  denotes  an  initial  homogeneous  shear  prestress,  R0  is  the  failure  zone  radius,  a  is  the  Poissons 
ratio  in  the  medium  and  p  the  rigidity.  Thus,  for  this  simple  case  for  which  wc  can  expect  to  obtain  the 
largest  gravity  induced  radiation  effect  associated  with  failure  in  a  stressed  medium; 
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At  a  frequency  of  1Hz  and  at  a  distance  of  about  one  wave  length  at  HIz  (ie.  abuui  5  kin.)  from  die  source 
origin,  we  have: 


C(G) 

Xo 


<  32  n2  p  G 


5(l-2q) 


7  -  5a 


Ro 


'pJ 


2p 


20 


(A-16.) 


Taking  a  large  failure  zone  such  that  R0  /  vp  =  1,  so  that  a  distance  of  one  wave  length  at  one  hertz  is  close 
to  the  failure  zone  boundary,  and  an  initial  shear  strain  such  that  |  x $/2ji  |  =  10~3,  which  is  typical  of  the 
level  producing  failure  ( ie .  hundred  of  bars  of  shear  stress  or  less),  then  with  a-  1/4  and  p~5  gm/cm3  as 
representative  values,  we  have 


i  i  <  lo-8 

where  this  value  of  |  Xo  I  has  the  units  and  the  magnitude  of  the  dynamic  strain  field,  in  this  case  near  the 
failure  zone  boundary.  Similarly,  the  spectral  density  of  the  rotation  components  in  (A-14b.)  will  be  such 
that 


1  x/G)  I  <10'8 

at  1  Hz  and  at  a  distance  of  one  wavelength. 

The  order  of  magnitude  bounds  for  these  quantities,  all  of  which  should  be  of  the  order  of  the  spec¬ 
tral  levels  of  the  dynamic  strain  changes  in  the  medium  if  they  are  lo  produce  the  seismic  radiation  from  a 
tectonic  source,  arc  clearly  many  orders  of  magnitude  less  than  those  observed.  In  particular,  we  can 
expect  them  to  be  of  about  the  order  of  (to'1)  |  t,(I3)/2  p  |  near  the  failure  zone,  which  m  the  present  case 
implies  a  level  of  about  10^*  at  1  Hz.  The  values  predicted  for  gravity  change  induced  effects  arc  there¬ 
fore  at  least  four  orders  of  magnitude  too  small,  as  would  be  expected.  The  same  analysis,  when  applied 
to  the  effects  of  gravity  changes  inside  the  failure  zone  alone  such  as  arc  implicit  in  the  "stress  glut"  for¬ 
mulation,  produce  similar,  but  even  smaller,  effects. 

The  origin  of  the  small  gravity  effect  is  associated  with  the  small  value  of  the  "seismo-gravity  cou¬ 
pling"  factor,  K.G  =  8rc2p2G,  appearing  in  the  basic  equivalent  force  field  relation  of  (A- 1 2. ).  Here, 
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G  =  6.67x10  8  (cgs  units)  and  so  KG  ~  KT4.  The  coupling  factor  therefore  scales  (down)  the  dynamic  field 
by  about  lO^1  and  produces  the  small  gravity  related  effects  just  described. 


By  way  of  contrast  and  in  order  to  show  the  role  of  the  coupling  more  precisely,  the  magnitude  of 
the  seismic  radiation  due  to  stress  relaxation  around  such  a  spherical  failure  zone  is,  to  the  same  approxi¬ 
mation  as  was  used  for  the  gravity  effect  calculation  (see  Archambeau,  1972  for  exact  results): 


1  (r  ,  w)  = 


=  5(1  -  2a) 

1  -5a 


—  hj2)(kpr)  P I  (cos  0)  cos  <f> 


with  similar  relations  for  the  rotation  field  components.  At  1  Hz  and  at  one  wave  length  from  the  source 
origin,  we  have  for  the  magnitude  of  the  dilatation  field: 


XoR)  (r ,  (0) 


5(1  -  2a)  Roj3  *13 
7  -  5o  j  [  vp  j  2p 


When  Rq/Vj,  is  near  unity  then  this  magnitude  estimate  ,  at  a  distance  near  the  failure  boundary  and  for 
a  =  1/4,  is  of  the  order  of  the  prestrain  given  by  |  x$/  2  p  | .  Thus,  the  dilatation  at  1  Hz  (and  likewise  the 
rotation  field)  is  of  the  order  of  the  prestrain,  or  of  the  order  1CT3,  if  we  use  the  same  prestrain  level 
assumed  in  the  gravity  effect  calculation.  We  also  observe  that  this  radiation  field  has  the  same  spatial 
dependence  and  scales  with  source  dimensions,  prestrain  and  medium  properties  in  the  same  way  as  the 
gravity  induced  effect,  although  its  frequency  dependence  is  somewhat  different.  This  is  as  would  be 
expected  since  it  L'  this  direct  radiation  effect,  associated  with  the  relaxation  of  the  prestress  field,  that 
produces  the  related  gravity  effect.  However,  the  ratio  of  the  magnitudes  of  the  two  effects  is  seen  to  be 
such  that: 


I  X6G)  I  /  I  Xo{R)  I  <  -y-  =  10-4 

at  1  Hz  and  at  all  spatial  locations  in  the  medium.  A  similar  relationship  holds  at  other  frequencies  as 
well,  as  can  be  verified.  However  a  demonstration  at  1  Hz  is  sufficient.  Thus,  the  gravity  effects  an' 


clearly  negligible  compared  to  direct  stress  relaxation  produced  radiation. 
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